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INTRODUCTION 


This  is  an  interim  progress  report  on  Contract 
Nonr  2798 (01 ) (FBM)  covering  the  work  done  on  submarine 
communications  at  the  University  of  New  Mexico  during 
the  last  three  months  of  I96I  and  the  calendar  year  1962. 

It  was  prepared  to  supplement  the  oral  presentation  to 
be  given  to  the  Navy  Department  at  the  Naval  Research 
Laboratory  on  February  I8  and  I9,  I965.  A  similar,  oral 
and  written,  report  was  rendered  to  the  U.  S.  Navy  on 
October  13/  1961.  It  summarized  the  progress  on  the 
contract  up  to  October,  I96I. 

Bimonthly  progress  reports  are  submitted  to  ONR 
except  for  the  months  of  July  and  August,  when  monthly 
progress  reports  are  rendered.  These  reports  give  a 
detailed  account  of  the  work  as  it  progresses.  The  purpose 
of  this  report  is  not  to  itemize  the  work  done  on  the 
project  during  the  period  reported  upon,  but  to  summarize 
and  highlight  the  work  on  the  tasks  on  which  more  signifi¬ 
cant  progress  was  made.  For  more  detailed  information 
reference  should  be  made  to  the  reports  published  on  the 
contract.  Such  list  appears  elsewhere  in  this  report. 

It  should  be  noted  that  this  is  a  progress  report  and, 
therefore,  the  results  and/or  conclusions  presented  should 
not  be  taken  as  final.  Changes  in  the  findings  may  have  to 
be  made  as  additional  information  and  data  become  available. 
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SUMMARY  OF  GENERAL  PROGRESS 


During  the  last  three  months  of  I96I  and  the  calendar 
year  1962  the  investigations  conducted  under  Contract 
Nonr  2798(01) (FBM)  continued  in  general  along  the  same 
lines  as  the  year  before.  However,  our  goals  became 
clearer,  our  tasks  better  oriented  and  our  studies  directed 
toward  the  solution  of  certain  more  defined  problems. 
Several  tasks  were  con^leted,  others  were  modified  and 
several  new  ones  initiated. 

From  the  general  approach  to  the  problem  of  investi- 
gating  submarine  communications,  the  work  at  UNM  on  the 
ELF/VLF  portion  of  the  frequency  spectrum  can  be  divided 
into  three  general  areas;  these  are: 

1.  study  of  ELF/VLF  propagation  in  air,  sea  and 
the  ionosphere. 

2.  Investigation  of  ELF/Vlf  antennas  in  these 
media . 

3.  Modeling  and  experimentation  to  confirm  results 
of  the  theoretical  investigations. 

Bach  of  these  areas  is  subdivided  into  a  number  of 
tasks  with  specific  goals.  A  general  summary  of  our  pro¬ 
gress  on  each  task  follows. 

1.  The  study  of  BLF/VLF  propagation  consisted  of  the 
following  tasks : 
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a.  Undersea  Propagation.  The  main  effort 
under  this  task  consisted  of  evaluating  the  effect  that 
a  rough  sea  surface  (wave  height  and  length)  has  on  the 
an^litude  and  phase  of  an  electromagnetic  wave  transmitted 
from  the  air  to  the  sea;  this  is  to  be  followed  by  a  study 
of  the  effect  such  a  surface  may  have  on  a  wave  going 
from  the  sea  to  the  air.  In  other  words,  the  purpose  of 
this  task  is  to  develop  a  technique  by  means  of  which  the 
fields  near  the  interface  can  be  calculated  when  the  fields 
some  distance  from  the  interface  are  known.  The  initial 
phase  of  this  task  consisted  of  defining  or  describing  the 
surface  of  the  sea,  including  its  various  states  of  roughness. 
Having  done  this,  the  second  phase  is  to  solve  the  boundary 
value  problem.  We  have  made  some  progress  on  this  task  and 

a  more  detailed  report  on  it  follows. 

b.  Subterranean  Propagation.  This  task,  whose 
goal  was  to  investigate  the  feasibility  of  long-range  low 
frequency  propagation  in  the  earth's  crust,  was  concluded. 
Based  on  experiments,  we  found  that  even  voider  the  best  of 
conditions  low  frequency  propagation  in  the  earth's  crust 
beyond  the  range  of  several  hvmdred  miles  was  not  feasible. 

We  will  give  a  concluding  report  on  this  task. 

c.  Satellite -Submarine  Propagation.  The  purpose 
of  this  task  is  to  investigate  ELF/Vlf  propagation  in  the 
ionosphere  in  order  to  determine  the  feasibility  of  propa¬ 
gating  ELF/VlF  signals  radiated  from  satellite  antennas. 
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The  initial  effort  consisted  of  studying  the  general  theory 
associated  with  propagation  in  anisotropic  media  such  as 
the  ionosphere.  Approximate  solutions  were  developed  for 
the  radiation  field  from  a  dipole  source  in  the  ionosphere. 
Then  power  relations  and  radiation  resistance  for  elementary 
dipoles  operating  at  VLF  in  the  ionosphere  were  derived. 

We  found  that  the  same  approach  applies  to  the  VHF  portion 
of  the  spectrum.  We  will  go  into  more  detail  on  the  progress 
on  this  task. 

d.  Study  and  Measurement  of  Noise  in  Sea. 

The  initial  effort  under  this  task  consists  of  finding  a 
suitable  method  to  evaluate  sensors  in  a  conducting  medium 
such  as  sea  water,  ^is  includes  the  coaucial  type  (both 
shorted  and  open),  the  loop  type  of  several  configurations 
and  the  toroidal  type.  Expressions  are  being  derived  for 
the  open  circuit  voltage  for  the  three  types  of  sensors. 

This  requires  an  expression  for  the  self-inqpedance  of  each 
type  of  sensor.  As  of  this  date  the  work  on  this  task  has 
not  progressed  sufficiently  to  report  on  it  in  more  detail. 

2.  Our  investigations  of  antennas  covered  five  separate 
tasks.  The  first  one  consisted  of: 

a.  Submerged  Anteimas.  Following  the  work  that 
Dr.  Anderson  did  several  years  ago  on  submerged  ELF  antennas, 
we  are  reconsidering  our  approach  to  this  task.  Anderson, 
in  his  study,  made  some  simplifying  assunptions  and  approxi¬ 
mations.  His  results  are  applicable  only  at  large  distances. 
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Anderson  felt  that  some  kind  of  a  numerical  process  would 
have  to  be  used  to  remove  these  limitations.  Further 
study  showed  that  the  whole  problem  was  much  more  compli¬ 
cated  than  at  first  anticipated.  It  appeared  that  integrals 
of  the  type  Anderson  used  occur  frequently  in  electromagnetic 
and  acoustic  problems  involving  layered  media.  Ohere  was  a 
need  to  evaluate  such  expressions  and  it  appeared  desirable 
to  find  a  computer  technique  to  make  the  solution  of  such 
integrals  routine.  After  some  study  of  the  problem,  the 
conclusion  was  reached  that  to  solve  these  integrals  for 
the  variation  of  the  several  parameters  involved  would  be 
an  enormous  task.  It  has  been  shown  that  without  the  iono¬ 
sphere  a  submerged  horizontal  electric  dipole  can  be  replaced 
by  another  horizontal  electric  dipole  at  the  air-sea  boundary. 
It  appears  that  such  an  approach  to  the  problem  may  be 
superior  and  more  direct  to  that  of  evaluating  the  integrals 
on  a  computer.  We  are  looking  into  this  approach  further. 

Under  this  task  we  have  also  continued  the  study  of 
coaxial  and  toroidal  receiving  antennas  in  sea  water.  This 
phase  of  the  task  will  be  covered  in  more  detail. 

b.  ELF  Antennas  for  Launching  Waves  in  Air.  A  sub¬ 
task,  which  consisted  of  the  study  of  the  feasibility  of 
using  an  island-sea  configuration  as  a  natural  slot  (aperture) 
antenna  for  operation  in  the  Elf  electromagnetic  spectrum, 
was  closed  out  in  I96I  when  it  became  apparent  that  others 
have  already  investigated  the  problem.  The  concluding  report 
on  this  phase  of  the  task  was  published  last  year  as  TM-6. 
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An  investigation  was  initiated  in  September  1962 
to  determine  if  certain  subsurface  antenna  configurations 
could  be  used  to  launch  ELF  waves  in  air  in  a  manner 
comparable  to  the  above -the -surface  vertical  antennas. 

The  configuration  now  under  consideration  is  a  toroid 
excited  by  one  or  more  current  windings.  This  task  has 
not  progressed  sufficiently  to  report  on  it  in  more  detail. 

c.  Air -to -Submarine  Antennas.  The  purpose  of 
this  task  is  to  evaluate  the  effectiveness  of  a  horizontal 
against  a  vertical  airborne  antenna  from  the  standpoint  of 
radiation  at  very  long  ranges.  The  study  has  centered  at 
15  kc/s.  The  problem  has  required  the  development  of  a 
technique  to  determine  field  strengths  at  long  ranges. 

AlsO/  the  effect  of  anisotropy  of  the  ionosphere  on  the 
propagation  of  VLP  waves  is  being  investigated.  We  will 
describe  our  progress  on  this  task  in  more  detail. 

d.  Satellite -to -Submarine  Antennas.  This  task 
is  the  second  phase  of  the  satellite-to-submarine  propaga¬ 
tion  study.  This  phase  has  centered  on  the  calculation  of 
approximate  expressions  for  the  fields  in  a  region  immediately 
adjacent  to  a  finite  antenna  structure  in  an  anisotropic 
homogeneous  ionosphere.  These  .will  be  used  in  an  attempt  to 
determine  the  driving -point  impedance  of  a  cylindrical  struc¬ 
ture  radiating  in  such  an  ionosphere.  We  will  go  into  more 
detail  on  the  progress  on  this  phase  of  the  task. 
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e.  VLF  Airborne  Antennas.  This  task  had  as  its 
goal  the  determination  of  the  advantages,  if  any,  that  a 
parallel  wire  side-loaded  transmission  line  may  have  if 
used  as  a  trailing -wire  antenna  and  to  suggest  a  design 
for  a  practical  form  of  this  type  of  an  antenna  for  use 
in  the  VLF  spectrum.  The  first  phase  of  this  task,  that 
is  the  evaluation  of  the  advantages  of  a  side -loaded 
trailing-wire  transmission  line  and  a  method  of  analysis 
of  such  an  antenna,  has  been  completed.  The  results  of 
the  experimental  work  on  a  scaled  model  were  quite  success¬ 
ful.  However,  the  second  phase,  the  design  of  a  practical 
and  flyable  antenna,  is  still  under  way.  We  will  present 
our  progress  on  this  task  in  more  detail. 

5.  The  third  area  of  our  work  was  in  experimentation, 
modeling  and  testing.  To  confirm  the  resul^ts  or  conclusions 
of  our  studies  some  experimental  data  is  required.  It  is  also 
needed  to  get  some  concept  of  how  the  dependent  parameters 
vary  when  the  independent  variables  are  changed.  Because  of 
the  tremendous  effort  that  would  be  involved  in  any  full-scale 
experimental  work  at  ELF/VLF,  any  work  of  this  type  that  we 
do  must  be  appropriately  scaled.  Our  work  of  this  type  is 
limited  in  scope  and  restricted  only  to  a  back-up  for  our 
investigations.  It  consists  of  laboratory  scaled  models  of 
antennas,  of  the  two-  and  multi-layer  problems  and  of  scaled 
in5>edance  measuring  circuitry.  We  will  cover  this  phase  of 
our  work  in  more  detail. 


THE  EFFECT  OF  THE  ROUGH  AIR-SEA  INTERFACE  ON 
ELF/VLF  ELECTROMAGNETIC  WAVE  PROPAGATION 

The  problem  of  E.M.  wave  propagation  near  the  rough  air- 
sea  interface  falls  naturally  into  two  parts: 

1.  The  description  of  the  air-sea  interface  (i.e.,  the 
sea  surface); 

2.  The  solution  of  the  boundary  value  problem  for  the 
E.M.  fields  in  the  presence  of  the  rough  interface. 

As  Part  One  is  necessary  in  order  to  define  clearly  the 
problem  considered  in  Part  Two,  the  present  state  of  knowledge 
about  the  mathematical  description  of  the  sea  surface  is  first 
briefly  reviewed.  Then  the  problem  of  Part  Two  is  considered 
in  light  of  this  description. 

The  Nature  of  the  Sea  Surface 

Even  casual  observation  shows  the  great  irregularity  of 
the  sea  surface;  no  single  wave  retains  its  identity  long; 
the  period,  form,  etc.,  vary  greatly  even  for  consecutive 
waves.  Indeed,  the  sea  surface,  in  a  rough  sea,  seems  to 
vary  almost  randomly  both  in  time  and  space. 

Presented  below  is  a  very  brief  review  of  the  present 
state  of  knowledge  about  the  nature  of  the  sea  surface.  How¬ 
ever,  the  study  of  the  sea  surface  continues  and  changes 
continue  to  be  made  in  the  theories  and  mathematical  formula¬ 
tions;  many  questions  about  the  sea  surface  still  remain  un¬ 
answered  and  no  known  theory  is  in  any  sense  complete. 
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The  study  of  the  sea  surface  can,  for  our  purpose,  be 
placed  into  three  general  categories: 

a.  The  study  of  the  sea  surface  by  classical  hydro¬ 
dynamics  . 

b.  The  study  of  the  sea  surface  by  statistical  methods. 

In  this  category  only  the  linearized  problem  will  be  considered 
as  this  allows  us  to  obtain  general  results.  By  assuming  a 
linearized  free -surface  boundary  condition,  the  problem  becomes 
linear  and  the  sea  surface  can  then  be  described  by  a  known 
random  process  (  stationary  Gaussian  process). 

c.  Other  theories  which  consider  the  non-linear  effects 
of  the  sea  surface.  In  this  case  the  results  obtained  are  far 
less  useful  for  a  description  of  the  sea  surface  than  the  re¬ 
sults  of  Category  b.  However,  such  results  place  the  limita¬ 
tions  of  Category  b  in  the  correct  perspective. 

Before  going  into  a  more  detailed  study  of  the  sea  surface 
a  few  general  comments  are  appropriate.  First,  the  present 
theories  assume  a  fixed  meteorological  condition.  Then  the 
description  of  the  sea  surface  given  by  the  theory  holds  only 
as  long  as  this  assumption  is  approximately  true.  A  theory 
based  on  the  correlation  of  the  changes  in  meteorological  con¬ 
ditions  with  changes  in  the  sea  surface  could  be  used  but  would 
be  somewhat  complex;  however,  in  some  cases  this  may  be  necessary. 
Second,  the  known  theories  consider  only  the  gross  meteorological 
conditions  (e.g.,  average  wind  velocity,  average  fetch,  etc.). 

In  practice  this  is  all  that  can  be  assumed  without  making  the 
problem  inordinately  complex.  Third,  no  current  theory  is 
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complete  in  the  description  of  the  sea  surface.  Under  certain 
meteorological  conditions  one  theory  may  be  approximately  cor¬ 
rect,*  however,  it  fails  when  meteorological  conditions  change. 

It  should  be  noted  that  no  theory  now  available  gives  the  com¬ 
plete  description  of  the  detailed  properties  of  the  sea  surface; 
only  the  gross  features  can  be  accurately  described  mathematic¬ 
ally.  A  brief  discussion  of  the  three  general  categories  follows 

a.  Classical  Hydrodynamics.  When  the  depth  of  the  sea  is 
large,  a  solution  of  the  hydrodynamic  equations  which  represent 
the  sea  surface  is  a  trochoidal  wave  [1,2].  The  parametric 
equation  of  the  trochoid  is: 

X  =  rO  -  a  sin  0 

z(x)  =  r  -  a  cos  G 

where 

r,  a  -  fixed  parameters  for  a  given  trochoid. 

0  -  parameter  that  varies  (i.e.,  generates  the  curve). 

The  trochoid  is  a  two-dimensional  wave,  which  "could  exist" 
in  a  swell.  Using  the  trochoid  to  describe  the  sea  surface 
makes  such  a  surface  regular  and  uniform  while,  in  fact,  the 
sea  surface  is  very  irregular  even  in  a  swell. 

A  more  general  model  of  the  sea  state  would  include  the 
irregularity  of  the  sea;  however,  such  a  formulation  is  too 
complex  for  practical  use.  We  consider  this  model  as  some 
of  the  results  will  be  used  later.  Leaving  out  the  details, 
the  following  equation  of  a  simple  harmonic  progressive  wave 
is  a  solution  to  the  linearized  hydrodynamic  equations: 
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5(x,y,t)  =  A  cos  (x  cos  0  +  y  sin  0)  -  t  +  ej  (l) 

where 

S  -  sea  surface 

A  -  amplitude  of  the  simple  harmonic  progressive  wave 

T  -  period  of  the  simple  harmonic  progressive  wave 

L  -  wavelength  of  the  simple  harmonic  progressive  wave 

c  -  velocity  of  the  simple  harmonic  progressive  wave 

0  -  direction  of  propagation  measured  with  respect  to 
the  +x  axis 

c  -  phase  atx=y=t=o  (arbitrary) 

If  we  let: 

p  =  XX  +  yy  -  position  vector  in  the  horizontal  plane 

Jc  *  k  X  +  k  y  -  propagation  vector  in  the  horizontal 
X  y 

plane 

2 

k  »  (lc|  *  -  wave  number  [J] 

y  ii 

g  -  acceleration  due  to  gravity 

0)  =  -  radian  frequency  of  the  wave 

then: 

5(p  »t)  =  A  cos  (ic’p  -  (ot  +  e)  .  (2) 

If  we  assume  that  the  waves  are  progressing  in  the  +x 

direction  (i.e.,  J'jj  >  0)  which  is  reasonable  if  the  wind  is  in 

the  +x  direction,  the  general  solution  is  [4]: 

Tr/2  f 

a(a),0)  cos  (Jc*p-a)t)  +  b(cD,0)  sin  (jc*p -6Dt)JdcD  d0 

(3) 

where 

a(cD,0)  and  b(a),0)  are  the  spectra  of  5(p  ,t) 


!(p  .t)  -J 


-V/2  « 
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If  S(0/y,t)  is  known,  we  can  obtain  a(cu,t)  and  b(a),t)  and 
therefore  $(x,y,t).  However,  it  is  clear  that  the  deterministic 
model  of  the  sea  surface  as  given  above  is  not  practical.  As 
has  been  the  case  with  many  problems  whose  conplexity  defy 
deterministic  solution,  one  next  atten^ts  to  formulate  the 
problem  in  probabilistic  terms.  We  will  now  consider  such  an 
atten^t. 


b.  Statistical  Description  of  the  Sea  Surface.  In  the 
past  few  years  there  has  been  an  increased  tendency  to  treat 
many  natural  phenomena  as  random  processes.  The  main  feature 
of  such  a  process  is  an  indeterminacy  in  the  expected  behavior 
of  a  single  occurrence  coupled  with  strong  statistical  proper¬ 
ties  for  a  large  number  of  occurrences.  The  "indeterminacy" 
in  the  sea  surface  comes  from  the  complexity  of  its  mathemat¬ 
ical  description  [e.g.,  finding  S(0,y,t)]. 

Chronologically,  experimental  data  first  lead  to  the 
assertion  that  the  sea  surface  could  be  represented  approxi¬ 
mately  as  a  stationary  multivariate  Gaussian  process.  From  a 
theoretical  view,  it  can  be  shown  by  using  equation  (2)  with 
the  assumption  that  the  random  variable  e(u>,0)  has  a  uniform 
distribution,  that  the  general  solution  for  the  sea  surface  is 


C4]: 

?(p  .t) 


JJ 

0  -V 


-  cot  + 


c(co,0)]  yA^(cD,0)da)  dO 


(4) 


This  expression  is  a  multivariate  Gaussian  process,  where 
A^(cu,0)  is  the  energy  density  spectrum.  This  ^as  does  (2)^ 
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assumes  linearized  equations  and  boundary  conditions;  if  these 
assumptions  hold,  the  surface  is  given  by  (4).  However,  as 
stated  above  the  experimental  data  implies  that  the  surface  is 
only  approximately  Gaussian,  the  error  being  due  in  large 
measure  to  the  nonlinear  effects  neglected  by  this  theory. 

The  same  results  can  be  given  in  a  standard  notation  such 
as  used  in  reference  [^] 

IT  00 

5(p*  /t)  =J  J'A(aj,0)  cos  (ic*p  -  ojt  +  e(aj,0))  dco  d0  (5) 

-IT  0 

where: 

1  2 

A  (£0,0)  =  e(£o,0)  is  the  energy  density  spectrum. 

Using  this  Gaussian  model  of  the  sea  surface,  many  of  the 
general  properties  of  the  "sea  state"  have  been  calculated. 
[4,3/6].  The  basic  question  to  be  resolved  is:  Can  the  multi¬ 
variate  Gaussian  process  describe  accurately  the  real  sea  sur¬ 
face?  If  not,  from  what  standpoint  is  it  deficient? 

To  answer  this  question,  we  consider  a  specific  model  of 
the  sea  surface  (i.e.,  a  given  spectrum  E((o, 0))  and  then  com¬ 
pare  this  random  process  with  the  known  properties  (experimental 
data)  of  the  sea  surface.  Naturally,  for  a  different  spectrum, 
the  properties  of  the  sea  predicted  by  the  random  process  will 
be  different.  The  best  known  spectrum  of  the  "sea  surface"  is 
a  semi-en^irical  expression  given  by  Neumann  for  a  fully 
developed  sea  [4]: 

E(ai)  to  -  ;  -g 


(6) 
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where  v  is  the  velocity  of  the  wind  "generating"  the  "sea". 
The  total  energy  E  for  this  spectrum  then  becomes: 

"  3/2 

E  =jE(a))  =  °T3/23  (7) 

0  ^ 


This  seems  to  be  in  agreement  with  some  experimental  data  if 
c  =  3.05  mVsec^.  To  study  in  detail  what  this  model  predicts, 
we  will  consider  the  shape  of  the  spectrum  and  some  of  its 
results . 

Neumann's  spectrum  rises  rapidly  at  oi  =  XT^v  ^ 

maximum  at  o)  =  ^  g/v.  For  large  cu,  e(co)  m  and  for  small 
2  2  2 

CO,  e(co)  w  e“^^  ^  ;  therefore,  there  is  very  little  energy 

in  the  low  frequency  (long  wavelength)  part  of  the  spectrum. 

Neumann's  spectrum  is  for  a  fully  developed  sea  only. 

In  actuality,  of  course,  the  sea  may  not  be  fully  developed. 

The  growth  of  the  sea  waves  depends  upon  the  fetch  (distance 
over  which  the  wind  blows),  the  duration  (length  of  time  the 
wind  has  been  blowing)  and  naturally,  the  wind  velocity. 

There  are  methods  that  consider  the  problem  of  duration -limited 
and  fetch-limited  seas  (but  only  approximately). 

Prom  en^jirical  data  an  approximate  directed  spectr\am 
(which  depends  on  the  direction  relative  to  the  wind  as  well 
as  on  the  duration  and  fetch  of  the  wave)  is  [4]: 


E((u,e)  =  < 


e 


-2gVa)^v^ 


cos  0  for: 


0) 
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where 

0  -  polar  angle  with  reference  to  the  x  axis  (the  wind  is 
is  blowing  in  the  +x  direction) . 

(0.  -  intersection  radian  frequency  (a  function  of  duration 
*  and  fetch) . 

The  above  model  is  at  best  an  approximation  to  the  sea 
surface  and  accurate  only  under  certain  conditions.  There  are 
better  models  [3,7];  they  are,  however,  more  con^lex  concerning 
the  directional  part  of  the  formula;  other  inaccuracies  have 
also  been  noted  in  all  the  models.  This  is  particularly  true 
for  high  sea  states  where  the  sea  may  not  be  fully  developed. 
Other  empirical  formulas  are  available  to  represent  the  sea 
surface  in  this  case  [8].  It  is  still  not  conclusively  known 
how  well  the  statistical  results  of  this  model  hold  for  actual 
sea  waves. 

Neumann's  model  of  the  sea  surface  contains  some  energy 
in  the  high  frequency  part  of  the  spectrum  and  predicts  the 
sea  to  be  "con?>letely  covered  by  short  wavelength  ripples." 
Physically,  this  seems  reasonable  even  though  the  above  spec¬ 
trum  is  known  to  be  deficient  in  the  high  frequency  range. 

The  error  in  the  high  frequency  part  of  the  spectrum  is  due 
to  the  fact  that  non-linear  effects  are  prevalent  at  these 
frequencies;  we  consider  this  effect  in  the  next  section. 

c.  The  Generation  of  the  Sea  Surface  and  Non-Linear 
Effects .  Phillips  C9]  has  described  the  generation  of  the 
sea  surface  by  considering  turbulent  fluctuation  of  the  air 
(wind)  above  the  air-sea  interface.  At  present  this  theory 
seems  to  be  accepted  as  correct  for  the  "original"  formation 
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of  the  sea  surface.  However,  once  the  waves  have  been  formed, 
their  growth  occurs  through  other  mechanisms. 

At  present  there  are  two  principal  theories  on  the  growth 
of  waves: 

(a)  Phillips'  resonance  theory,  [10] 

(b)  Miles'  sheer-flow  Instability  theory,  [11] 

These  two  theories  do  not  agree  In  many  respects  (e.g., 

Phillips'  theory  gives  a  "growth  rate"  proportional  to  time 
and  Miles  gives  a  "growth  rate"  exponential  with  time).  While 
both  mechanisms  play  a  role  In  the  growth  of  waves,  Phillips 
has  recently  given  a  description  of  the  domains  of  dominance 
of  each  mechanism.  These  theories,  while  they  do  give  some 
understanding  of  the  sea  surface,  are  unfortunately  of  no 
qualitative  help  at  present  In  describing  the  sea  surface. 

We  will  consider  only  one  other  mechanism  of  energy  trans¬ 
fer,  the  breaking  of  the  waves.  Under  eufflclently  high  wi^ds, 
the  energy  transfer  from  the  breaking  of  waves  reaches  an 
equilibrium;  then  for  that  range  of  frequencies  where  the 
non-llnear  (breaking)  effects  are  Important  the  energy  spectrum 
Is  given  by: 

E(a))  ■  ag^o)"^  where  a  ■  10”^  (empirical  constant) 

The  same  basic  conclusions  (under  basically  the  same  restrictions) 
were  obtained  by  Mikhailov  [12],  who  used  the  theory  of  turbu¬ 
lence.  There  Is  close  agreement  between  these  theories 
[Phillips  has  E((o)  w  co~^,  Mikhailov  cites  experimental  work 
which  gives  E(a))  m  to  a>~^;  both  theories  are  within  these 
bounds.  Note  Neumann  in  this  range  has  E((i>)  m  cd~^]. 
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Tick  [15,14]  has  considered  the  problem  of  non-linear 
effects  from  a  perturbation  point  of  view.  He  adds  to  the 
"linear  (Gaussian)  spectrum"  a  correction  due  to  the  non¬ 
linear  effects.  However,  the  statistics  of  the  "non-linear" 
spectrum  are  not  known. 

All  theories  (as  those  above)  which  do  not  directly  con¬ 
sider  capillary  waves  (waves  of  very  short  wavelength)  do  not 
hold  for  frequencies  in  the  "capillary"  range.  While  much 
work  on  capillaries  has  been  done,  at  present  there  is  nothing 
available  on  correlating  this  theory  with  the  overall  sea  sur¬ 
face.  And  there  is  no  correct  information  available  on  the 
'bapillary" range  of  the  energy  density  spectrum.  Work  on  radar 
return  and  light  reflection  [15 ]  from  the  sea  surface  has  shown 
that  the  slopes  of  the  sea  surface  are  nearly  Gaussian  and  that 
the  sea  surface  curvature  is  highly  non -Gaussian.  The  usual 
experimental  data  taken  for  description  of  the  sea  surface 
does  not  include  frequencies  in  the  "capillary"  range. 

The  complexity  of  the  air-sea  interface  alone  is  enough 
to  make  one  resort  to  statistical  analysis;  however,  the  fact 
that  both  the  linear  model  and  the  experimental  data  give  a 
Gaussian  process  for  the  gross  features  of  the  sea  surface  im¬ 
plies  that  a  useful  model  would  be  a  statistical  one. 

If  one  is  willing  to  neglect  the  non-linear  effects  of 
the  sea  surface  and  use  the  Gaussian  model,  the  Ergodic 
theorem  [4]  iitplies  that  the  ensemble  analysis  of  statistical 
properties  is  the  same  as  space  (or  time)  analysis  of  the 
statistical  properties  (note  that  this  requires  a  stationary 
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process,  i.e.,  an  equilibrium  state).  There  are  questions  as 
to  the  validity  of  the  Gaussian  model  as  regards  the  statistical 
properties  of  the  sea  surface.  However,  if  one  is  only  interested 
in  the  gross  statistical  features,  then  the  Gaussian  description 
may  be  adequate.  However,  such  description  will  give  only  the 
"smooth  average"  shape  of  the  sea  without  any  details  such  as 
ripples  or  minor  variations.  The  experimental  data  on  sea  sur¬ 
face  variation  is  meagre;  the  processing  of  such  data  for  mean¬ 
ingful  results  is  long  and  costly;  therefore,  it  is  doubtful 
if  an  accurate  description  of  the  sea  surface  can  be  readily 
obtained.  The  shape  of  the  surface  of  the  sea  depends  upon 
past  as  well  as  upon  present  conditions;  however,  to  a  degree 
of  approximation,  past  effects  though  still  significant  may  be 
discounted.  Also  by  describing  the  sea  condition  by  means  of 
only  a  few  average  parameters,  the  sea  surface  will  not  be 
represented  accurately  in  every  detail,  but  from  the  practical 
approach,  this  may  be  all  that  can  be  done.  This  is  the  point 
of  departure  for  the  study  of  propagation  of  electromagnetic 
waves  near  the  air-sea  interface. 


The  Boundary  Value  Problem  for  E.M, 

Fields  in  the  Presence  of  the  Rough  Air-Sea  Interface 


"In  the  mathematical  sense  the  problem  at  hand  is  extremely 
complex,  since  it  is  impossible  to  use  the  method  of  separation 
of  variables  to  obtain  a  solution  of  the  wave  equation  which 
satisfies  the  boundary  conditions  specified  on  the  uneven  sur¬ 
face.  The  methods  for  obtaining  an  accurate  solution  of  this 
problem  in  general  form  has  (sic)  not  yet  been  found.  However, 
there  are  a  number  of  theoretical  papers  (which  have  appeared 
essentially  during  the  last  3-6  years)  in  which  a  number  of 
approximate  methods  for  confuting  the  field  have  been  developed." 
This  is  a  quotation  from  a  survey  paper:  "Theory  of  the  Scatter¬ 
ing  of  Waves  at  Periodically  Uneven  Surfaces,"  by  lu.  P.  Lysanov 
Cl6].  This  paper  treats  for  the  most  part  the  problem  of  a  scalar 
field  in  the  presence  of  a  periodically  uneven  (rough)  surface 
with  "natural"  boundary  conditions.  When  the  boundary  conditions 
of  the  actual  problem  are  not  the  "natural"  ones,  the  "natural" 
boundary  conditions  are  used  as  an  approximation.  The  "interface" 
boundary  conditions  are  not  considered.  The  problem  treated  is 
then  the  simplest  problem  considering  a  rough  surface;  even 
then,  this  problem  is  not  solvable  using  known  techniques.  Al¬ 
though  work  on  this  problem  has  continued  for  a  number  of  years, 
little  advancement  has  been  made;  no  new  methods  or  theories 
have  resulted;  and  the  investigation  has  generally  been  confined 
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to  refinement  of  techniques  given  in  Lysanov's  paper  and 
experimental  data.  The  theory  (techniques  of  "solution")  for 
the  rough  surface  {or  interface)  problem  is  reviewed  with  the 
discussion  limited  to  the  inadequacy  of  the  "solutions". 

Basically  Lysanov's  approach  will  be  followed.  He  con¬ 
siders  six  general  classifications  (or  techniques): 

1.  Rayleigh's  Method.  This  is  the  oldest  approach  and 

due  to  Rayleigh.  Rayleigh  assumed  that  the  reflected  (scattered) 
field  could  be  represented  as  a  sum  of  outward  (from  the  surface) 
directed  plane  waves.  However,  this  is  only  true  for  the  fields 
above  the  highest  point  of  the  surface.  Lippmann  Cl?]  showed 
that  the  assumption  of  Rayleight  was  not  quite  correct; 
unfortunately,  he  could  not  obtain  aun  "accurate"  solution  to 
the  problem  by  the  use  of  his  variational  technique. 

Barantsev  [ l8]  has  also  obtained  the  same  result  by  use  of 
Laplace  transforms.  In  Barauitsev's  approach  to  the  problem, 
as  in  Rayleigh's,  aui  infinite  number  of  algebraic  equations  must 
be  solved.  In  the  case  where  the  irregularities  are  small  (the 
maximum  height  of  the  surface  wave  is  much  smaller  than  a  wave¬ 
length  of  the  radiation  field)  Rayleigh  obtained  a  solution  for 
the  first  mode  (specular  reflection);  this  agrees  quite  well 
with  experimental  data  in  the  "far  field",  however,  near  the 
Interface  many  evanescent  (surface)  waves  are  present  and  an 
accurate  solution  is  not  practicable  (even  if  Rayleigh's 
assxomption  were  to  be  considered  correct)  [l9»25]. 

2.  The  Method  of  Small  Perturbations.  The  boundary  condi¬ 
tions  are  specified  on  the  uneven  surface,  z  =  5(x),  and  then 
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are  transferred  to  the  plane  z  =  0  by  means  of  an  expansion  into 
a  power  series  with  respect  to  5  •  It  has  been  shown  by  Lysanov 
that  this  method  leads  to  solutions  identical  to  Rayleigh's 
results.  It  was  used  by  Feinberg  for  solving  the  problem  of 
propagation  of  radio  waves  along  the  earth's  surface  and  more 
recently  by  Senior  [20]  who  considered  the  effect  of  surface 
roughness  on  the  reflection  of  E.  M.  waves.  Senior  was  able 
to  replace  surface  roughness  by  a  change  in  the  impedance  at 
the  boundary.  The  surface  impedance  is  then  a  function  of 
position.  Hessel  [21]  and  others  have  also  considered  the 
problem  of  varying  surface  impedance;  Hessel  seems  to  account 
for  Wood's  anomaly  this  way;  however,  this  method  is  useful 
only  in  the  far  field.  Bass  [22]  also  uses  the  method  of 
perturbations  to  obtain  results  parallel  to  Senior,  these  again 
are  valid  only  in  the  far  field.  Lysanov  [23]  ,  has  extended 
this  theory  to  point  sources. 

3.  The  Method  of  L.  M.  Brekhovskikh .  This  is  the  method 
of  physical  optics  [24]  .  As  in  this  method  the  fields  on  the 
surface  are  assumed,  the  problem  becomes  one  of  evaluation  of 
an  integral  [25]  which  cannot  be  readily  accomplished  for  the 
region  near  the  surface. 

A  related  approach,  the  Kline -Luneberg  theory,  has  been 
recently  used  to  calculate  the  fields  on  a  rough  dielectric 
surface  [26]  .  A  special  case  of  the  Kline-Luneberg  theory, 
referred  to  as  geometric  optics,  has  been  applied  by  Leontovich 
to  the  case  of  a  highly  conducting  medium  in  which  the  wave¬ 
length  is  "small,"  and  then  extrapolated  to  the  problem  of  a 
rough  interface  [27]  • 
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For  the  frequency  range  that  we  are  interested  in  (ELFA’LF) 
the  use  of  geometric  or  physical  optics  does  not  appear  generally 
valid. 

4.  The  Integral  Equation  Method.  The  problem  of  the 
scalar  boundary  value  problem  with  "natural"  boundary  conditions 
can,  by  use  of  Green's  theorem,  be  formulated  conveniently  as 
an  integral  equation.  This  equation  is  exact;  however,  it 
cannot  be  solved  without  making  some  approximations.  If  the 
rough  surface  satisfies  certain  conditions  (is  not  too  rough), 
an  approximate  integral  equation  is  obtained,  which  can  be 
solved.  There  is  a  difference  in  the  field  assumed  on  the  s'ur- 
face  by  geometric  optics  and  the  field  on  the  surface  obtained 
by  the  integral  equation  method  [28].  Unfortxinately,  the  solu¬ 
tion  to  the  integral  equation  is  not  at  present  known  to  be 
accurate  enough  for  our  purpose. 

5«  The  Method  of  Images.  The  method  of  images  can  be 
used  to  investigate  the  fields  in  the  presence  of  an  uneven 
surface  with  a  sufficiently  single  shape;  this  has  been  done  ly 
!^ersky  for  a  perfectly  reflecting  plane  covered  with  half¬ 
cylinders  or  hemispheres  (with  little  interaction  between 
scatterers).  Biot  [29]  considered  the  perfect  conduction  pis'loe 
covered  by  hemispheres  with  strong  interaction  (as  there  would 
be  in  tlie  case  of  the  sea  surface).  Twersky  also  has  considerr^d 
multiple  scattering  in  a  very  general  way.  Hiis  method  could  le 
used  to  find  the  fields  on  the  surface,  assumed  to  be  a  plar^^^ 
covered  with  hemispheres;  however,  to  obtain  a  more  realistic 
model  of  the  sea  surface  would  involve  higher  multi -pole  evp^-,;:.- 
sions  and  this  does  not  seem  practical  (though  Biot  thought  t:: 
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method  could  be  used  for  the  air-sea  interface  problem).  It 
does  give  a  good  idea  of  what  is  happening  to  the  E.M.  field 
in  the  air  near  the  air-sea  interface  (where  Biot  uses  a 
static  solution). 

6.  Method  of  Matching  Fields.  This  method  can  be  used 
only  in  the  case  where  the  rough  surface  is  such  that  the  space 
may  be  separated  into  regions  in  which  the  wave  equation  allows 
solution  by  the  method  of  separation  of  variables  (when  the 
wave  equation  is  written  in  an  appropriate  coordinate  system). 

As  this  is  not  possible  for  the  air-sea  interface,  we  will  not 
consider  this  method  further. 

Lysanov  references  several  papers  on  experimental  work, 
since  then  papers  covering  similar  work  have  been  piiblished 
which  seem  to  imply  Rayleigh's  theory  (including  the  second 
mode)  is  approximately  correct  for  the  low  frequency  problem 
except  near  =  Ay,  (A^,  wavelength  of  radiation;  Ay,,  wavelength 
of  the  surface)  C19]  .  Since  1958  many  papers  on  the  statistical 
analysis  of  the  reflection  of  sound  from  the  rough  sea  surface 
have  been  published;  these  use  either  Rayleigh's  assumption  or 
Kirchhoff's  approximation  for  the  solution  of  the  boundary  value 
problem  and  as  such  are  only  useful  in  the  far  field  [50,51]  . 

One  of  the  few  papers  considering  "interface"  boundary  con¬ 
ditions  on  an  irregular  surface  was  published  by  Walt  [52]  in 

» 

1959 >  this  paper  has  been  used  to  calculate  the  effect  of  the 
rough  surface  on  the  E.M.  fields  in  the  sea.  Walt  bases  his 
approach  on  the  work  of  Leontovich  [27]  and  obtains  the  field 
on  a  plane  in  the  sea;  he  then  solves  the  wave  equation  in  the 
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sea  with  this  as  a  boundary  condition.  A  short  discussion  of 
this  work  is  given  in  a  paper  by  Whalen  [52]  .  It  is  noted  there 
that  the  model  of  the  sea  as  used  in  Wait's  paper  was  not  too 
realistic.  Because  of  his  approximation.  Wait  uses  a  "local 
boundary  condition"  which  is  not  exact. 

Another  paper  of  interest  is  by  Bass  [5^]  >  considers 
non-local  boundary  conditions,  but  uses  a  perturbation  method 
which  invalidates  the  theory  near  the  interface.  He  also 
considers  the  effect  of  the  variation  with  time  of  the  rough 
surface . 

Lastly,  we  note  a  number  of  papers  on  radar  scattering 
from  rough  surface  including  the  surface  of  the  sea  as  given 
by  Lysanov.  These  usually  use  Kirchhoff's  approximation  which 
is  reasonable  in  this  frequency  range.  An  excellent  discussion 
of  the  radar  problem  is  given  by  Kerr  [2]  .  The  most  frequently 
used  method  is:  "The  construction  of  the  phenomenological 
theory  starting  from  the  statistical  properties  of  the  received 
signal."  [55] 

The  basic  limitation  of  currently  known  theories  of  propa¬ 
gation  of  E.M.  waves  in  the  presence  of  a  rough  interface  is 
that  no  theory  holds  accurately  near  the  interface.  Walt's 
approach  applies  the  "non  local"  effects  after  the  waves  are 
below  the  interface  and  not  at  the  interface;  this  introduces 
some  inaccuracy.  The  use  of  Kirchhoff's  approximation  at  ELF 
and  VLF  can  be  questioned.  No  known  theory  can  be  considered 
sufficiently  accurate  at  the  interface,  so  that  the  fields  on 
the  svurface  can  be  used  as  sources  Cfollowlng  Huygen's  Principle) 
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for  the  fields  in  the  sea.  Even  if  we  knew  the  fields  on  the 
interface,  the  use  of  Huygen's  Principle  in  its  usual  form 
may  not  be  very  practical  (could  involve  a  large  amount  of 
numerical  integration). 

As  noted,  the  known  attempts  to  solve  the  problem  for  the 
field  near  a  rough  fixed  surface  in  analytic  form  have  not  been 
very  successful.  The  question,  if  any  of  the  many  approximate 
methods  (of  which  the  method  of  small  perturbation  seems  most 
popular  at  present)  are  "sufficiently  accurate"  near  the  surface, 
cannot  be  answered  categorically  as  no  bounds  can  be  set  on 
the  errors.  Efforts  to  correct  this  have  not  been  successful 
and  it  is  felt  that  little,  if  any,  success  would  be  made  by 
continuing  this  approach.  While  further  work  on  a  general 
analytic  solution  to  this  problem  has  continued,  no  progress 
can  be  reported.  In  light  of  the  quotation  from  Lysanov,  it 
seems  that  other  methods  of  solution  should  be  attempted 
collaterally  with  a  general  analytic  study.  In  the  past, 
numerical  solutions  to  a  nximber  of  special  problems  have  furn¬ 
ished  material  for  a  general  theory.  In  this  case,  the  numerical 
solutions  to  some  special  cases  may  give  enough  insight  to  do 
the  seune  (i.e.,  give  intuitive  or  the  "needed  accurate 
assiimption" ) .  This  would  also  make  it  possible  to  compare  the 
results  from  the  many  approximate  methods  with  the  "correct" 
result  and  obtain  bounds  on  the  error.  Numerical  solutions 
could  also  give  the  "general  effect"  of  surface  roughness. 

It,  therefore,  seems  reasonable  to  turn  to  other  methods 
to  solve  the  problem  numerically.  A  short  outline  of  such 
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methods  is  available  in  literature  [36]  .  The  use  of  the 
following  two  general  methods  has  been  considered: 

1.  Analog  (simulation) 

2.  Numerical  (finite  difference) 

An  extensive  study  of  analog  simulation  is  given  by 
Karplus  [373  •  The  use  of  analog  computers  does  not  seem 
suited  to  the  problem.  However,  the  use  of  network  analogs, 
based  on  the  finite  difference  approximations,  does  appear  to 
apply.  Because  of  its  cost,  this  method  will  be  considered 
only  after  other  methods  are  found  unapplicable. 

The  numerical  solution  of  differential  equations* has  been 
under  extensive  study  for  the  last  several  years.  This  has 
produced  many  excellent  papers  and  books  [38,39]  on  the  subject. 
There  are  many  techniques  for  rapidly  obtaining  solutions  to 
boundary  value  problems;  however,  these  methods  are  not 
applicable  near  a  rough  interface.  At  present  our  calculations 
are  limited  to  two-dimensional  problems  (the  solution  for  one 
scalar  function).  Even  in  this  case,  the  work  involved  is 
extremely  long. 


CONCLUSION 

While  the  available  mathematical  descriptions  of  the  sea 
surface  are  not  accurate  at  very  long  and  very  short  lengths  of 
sea  waves,  the  available  expressions  for  energy  density  spectra 
give  the  correct  gross  features  of  the  sea  surface.  When  using 
the  finite  difference  method  for  the  determination  of  the 
electromagnetic  fields  near  a  rough  sea  surface,  only  the  gross 
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features  of  the  sea  surface  are  required.  It  appears  “best"  to 
use  "average"  values.  The  finite  difference  method  appears 
at  present  to  be  the  best  mathematical  approach  to  the  rough 
interface  problem.  Use  of  this  method  does  not  require 
detailed  knowledge  of  the  very  short  wavelength  ripples  or 
capillary  waves.  It  also  appears  reasonable  to  assume  that  the 
statistical  properties  of  the  gross  features  of  the  sea  surface 
are  given  by  a  joint  Gaussian  distribution  for  low  sea  states; 
for  high  sea  states  the  non-linear  effects  predominate  and 
the  statistics  are  not  known. 

The  first  case  considered  using  th©  finite  difference 
method  was  a  two-dimensional  sinusoidal  surface.  The  result 
agrees  with  other  theoretical  studies  [32]  near  the  surface 
or  the  water.  Using  this  method,  calculations  will  be  made 
to  see  if  the  seune  conclusion  applies  to  rougher  surface  under 
more  general  assumptions.  The  method  can  in  theory  be  used 
for  three  dimensional  problems;  we  will  see  if  this  can  be  done 
practically.  Such  an  approach  to  the  problem  will  require  some 
correlation  with  experimentally  determined  data. 
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Some  Useful  Formulae 

For  reference  some  useful  relations  pertaining  to  the  sea 
surface  are  given  below. 

If  S(r)  =  A  cos  (Jc*r  axt)  Is  a  solution  to  the  linearized 

hydrodynamic  equations,  then: 

2 

0)  ■  gk  tanh  kh 
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We  then  have  for  a  single  harmonic  wave  In  deep  water: 
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For  relatively  "regular" seas  (after  they  have  consolidated 
Into  a  "regular"  series  of  connected  troughs  and  crests). 
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ELECTROMAGNETIC  RADIATION  IN  MAGNETO-IONIC 
MEDIA  FOR  VLP  AND  VHF  SPECTRA 

A  method  of  analysis  is  presented  for  obtaining  approxi¬ 
mate  solutions  to  Maxwell's  equations  pertaining  to  current 
sources  in  an  homogeneous  magneto-ionic  medium.  The  analysis 
is  applied  to  the  radiation  of  electromagnetic  energy  in  the 
VLP  and  VHP  spectra  Cl]. 

The  analysis  is  directed  toward  the  development  of  a  set 
of  potential  functions  for  describing  the  electromagnetic 
fields  of  finite  electric  current  sources.  The  gauge  condi¬ 
tions  associated  with  the  potential  functions  are  discussed 
and  the  limitations  of  a  particular  gauge  condition,  the  tensor 
Lorentz  gauge,  are  discussed  in  detail.  The  analysis  is  then 
restricted  to  specific  frequency  spectra  and  regions  of  the 
ionosphere  which  are  applicable  for  communications  to  and  from 
space  vehicles  within  such  regions .  It  is  found  that  approxi¬ 
mate  solutions  of  Maxwell's  equations  may  be  obtained  from  a 
potential  function  which  satisfies  the  tensor  form  of  Helmholtz's 
equation.  Such  solutions  are  applicable  to  current  sources  in 
the  VLP  and  VHP  spectra  in  or  above  the  F-layer  of  the  iono¬ 
sphere  (i.e.,  equal  to  or  greater  than  1^0  kilometers  above 
the  earth's  surface) . 

Closed  form  solutions  for  the  potential  function  are  then 
derived  for  both  electric  and  magnetic  dipoles  in  a  homogeneous 
ionosphere  by  use  of  Fourier  analysis  cuid  the  associated  Green's 


functions . 
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Results  are  presented  for  the  radiation  characteristics 
of  electric  and  magnetic  dipoles  arbitrarily  oriented  with 
respect  to  the  static  magnetic  field  or  the  gyrotroplc  axis 
of  the  magneto- ionic  medium.  Main  emphasis  is  placed  on  the 
field  structure  in  the  radiation  zone.  From  the  field  struc¬ 
ture  the  Poynting  vector  is  determined  for  the  VLP  and  VHF 
spectra,  and  hence  the  net  radiated  power  from  the  dipole. 
Radiation  phenomena  are  compared  with  those  for  free  space 
for  the  same  dipole. 

The  first  general  solution  for  the  radiation  of  a  finite 
current  distribution  in  a  hcxnogeneous  anisotropic  medium  was 
given  in  an  elegant  presentation  by  F .  V.  Bunkin  [2].  Bunkin 
considered  only  monochromatic  source  functions  and  assumed 
the  anisotropic  medium  was  characterized  by  a  general  permit¬ 
tivity  tensor,  the  permeability  being  a  known  scalar.  He 
then  considered  only  the  partial  differential  equation 
satisfied  by  the  electric  field  vector  as  required  by  Maxwell's 
equations: 

7(7.?)  -  7^?  -  k^€?  =  ?  (1) 

e:  complex  dielectric  tensor 

k:  free  space  propagation  constant  in  CGS  units, 

specified  current  distribution 
c:  speed  of  light 

Utilizing  the  linearity  of  the  above  differential  equation, 
Bunkin  assumed  that  the  electric  field  vector  could  be  ex¬ 
pressed  in  the  Integral  form: 
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t  =  avj  (2) 

''i 

where  ?(r,rj^)  is  a  tensor  of  second  rank  depending  on  the  coor¬ 
dinates  of  the  point  of  observation  (r)  and  the  source  point 

(?l). 

In  a  concise  but  rigorous  manner  he  then  achieved  the 
formal  solution  for  the  electric  field  of  a  general  current 
distribution  in  a  homogeneous  anisotropic  medium.  He  then 
directed  his  attention  to  a  magneto-ionic  medium  such  as  the 
ionosphere  for  evaluation  of  the  integrals. 

Restrictions  of  Bunkln's  Solution 

Because  of  the  complexity  of  the  integral  solutions, 

Bunkin  was  forced  to  restrict  his  attention  to  the  radiation 
zone,  where  asymptotic  methods  could  be  applied  to  obtain  a 
formal  evaluation  of  the  integrals.  This  procedure  produced 
a  multipole  expansion  of  the  electric  field  due  to  an  arbi¬ 
trary  current  distribution  in  a  magneto-ionic  medium. 

Although  his  formal  solution  was  valid  for  euiy  current 
distribution  and  any  point  of  observation,  numerical  evalua¬ 
tion  of  such  a  solution  was  restricted  to  the  electric  dipole 
distribution.  Furthermore,  expressions  for  the  electric 
field  vector  were  given  in  terms  of  a  parameter  Introduced 
through  the  use  of  the  method  of  steepest  descent  (the  saddle 
point) ,  whose  value  is  not  easily  determined  for  any  medium 
and  point  of  observation.  More  specifically,  the  parameter 
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is  a  function  of  the  angle  between  the  observation  point  and 
the  static  magnetic  field  as  well  as  the  elements  of  the 
permittivity  tensor. 

In  conclusion,  Bunkln  considered  only  two  orientations 
of  the  dipolej  and  for  each  orientation  he  discussed  only  the 
solutions  at  particular  observation  points  which  produced 
simplified  values  of  the  parameters  Introduced  in  his  analysis. 

Although  Bunkin' s  presentation  was  quite  elegant,  it  did 
not  manifest  itself  as  a  working  paper  for  further  use  in  the 
radiation  characteristics  of  an  electric  dipole  in  an  aniso¬ 
tropic  medium. 

During  the  year  196O  Kogelnik  and  Kuehl  published  certain 
results  for  the  radiation  characteristics  of  an  electric  dipole 
in  an  anisotropic  medium  which  utilized  the  formulation  of 
Bunkin  [3,4].  Kogelnik  determined  the  quadratic  form  for  the 
radiated  power  of  an  electric  dipole  with  specific  orienta¬ 
tions  relative  to  the  static  magnetic  field.  Through  the  aid 
of  electronic  computers,  he  was  able  to  obtain  numerical  data 
for  the  power  relations  of  cui  electric  dipole  but  did  not 
explicitly  evaluate  the  electromagnetic  field  structure. 

Kuehl 's  paper  was  essentially  a  reworking  of  Bunkin 's  pre¬ 
sentation  with  the  exception  that  it  was  developed  in  much 
greater  detail  and  did  provide  asymptotic  expressions  for  the 
far-zone  fields  of  an  electric  dipole  in  an  homogeneous  aniso¬ 
tropic  plasma.  However,  Kuehl 's  results  were  also  restricted 
to  two  particular  problems. 
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His  first  result  was  a  highly  theoretical  situation  in 
which  it  was  assumed  that  the  external  magnetic  field  was 
infinitely  large  such  that  the  integral  solutions  were  con¬ 
siderably  simplified.  This  condition  does  not  occur  in  the 
ionosphere  and  hence  is  not  applicable  to  any  physically 
realizable  problem  of  radiation  in  the  earth's  atmosphere. 

Kuehl's  other  result  was  for  a  magneto-ionic  medium  in 
which  the  normalized  plasma,  normalized  gyro  and  normalized 
collision  frequencies  were  much  less  than  one.  Physically 
this  corresponds  to  source  functions  operating  in  the  VHP 
spectrum  and  distauices  above  the  earth's  surface  exceeding 
150  kilometers .  Indeed  this  result  is  applicable  to  com¬ 
munication  problems  for  space  vehicles . 

The  Need  for  Further  Investigation  of  the  Radiation  Problem 

Recent  experimental  work  in  the  ionosphere  has  required 
a  knowledge  of  radiation  phenomena  in  the  VLF  spectrum. 
However,  a  review  of  the  literature  reveals  that  little  in¬ 
formation  is  available  on  the  radiation  of  electromagnetic 
energy  in  the  VLF  spectrum  as  applied  to  the  ionosphere .* 
Theoretical  investigations  have  essentially  been  restricted 
to  those  discussed  above  which  are  not  readily  applicable  to 
the  VLF  spectrum.  A  detailed  investigation  of  the  work  of 
Bunkin  and  Kuehl  revealed  that  the  application  of  this  type 

* 

The  paper  by  R.  Mlttra,  "Solution  of  Maxwell's  Equations 
in  a  Magneto-Ionic  Medium  with  Sources,"  U.  of  Ill.,  had  not 
yet  become  available  when  K.  R.  Cook  finished  his  work. 
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of  analysis  to  the  VLF  case  is  prohibitive  in  view  of  the 
complex  integrations  required.  In  view  of  these  conclusions, 
and  the  need  for  information  pertaining  to  radiation  phenom¬ 
ena  in  the  VLP  spectrum,  it  was  apparent  that  a  different 
method  of  analysis  was  required  if  useful  results  were  to 
be  obtained. 

The  following  approach  describes  a  new  method  of  analysis 
for  obtaining  approximate  solutions  for  both  the  electromag¬ 
netic  fields  eind  radiated  power  for  current  sources  in  the 
VLP  and  VHP  spectra  in  a  magneto-ionic  medium. 

Considering  only  monochromatic  source  functions,  Meucwell's 
equations  for  anisotropic  media  may  be  expressed  by  the  fol- 


lowing  system: 

V  X  +  lOJIX^S  =  ^ 

(3) 

7  X  it  -  iu3€_K^  =  d 

o 

(4) 

5  =  tjit,  t  = 

(5) 

v-S  =  p,  v-S  =  0  . 

(6) 

Thus  we  see,  as  in  isotropic  media,  the  magnetic  field  inten¬ 
sity  represents  a  solenoidal  vector  field.  Therefore,  we 
may  conclude  that  there  exists  a  vector  ^  such  that 

^  =  V  X  X  .  (7) 

Prom  Maxwell's  equations  it  follows  that  the  electric  field 
is  representable  as, 

t  =  -  (icuix^^  +  vcp)  . 


(8) 
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The  potentials  X  and  cp  must  satisfy  the  differential 
equations, 

=  -  ?  +  77.x  4-  Icoe^icTcp  (9) 

7.(lC7Cp)  +  10511^7  .(  ic a)  =  -  p/e^  (10) 

We  now  assume  that  X  and  cp  satisfy  the  following  gauge  condi¬ 
tion, 

7(7-X)  +  l(J0e^l?7Cp  =  ^  (11) 

However,  we  find  the  assumed  gauge  condition  Imposes  restric¬ 
tions  on  the  vector  potential  X.  To  show  this  we  premultiply 
(11)  by  the  Inverse  permittivity  matrix  S,  to  produce  the 
result, 

Jt7(7.X)  (12) 

Therefore,  It  follows  that  the  vector  potential  must  satisfy 
the  condition, 

7  X  [5|7(7.X)  ]  =  ^  (13) 

From  the  above  equations  we  find  that  under  the  assumed  gauge 
condition  the  vector  potential  must  also  satisfy  the  equation 

7^X  +  kJi?X  =  -  J  (14) 

o 

In  general  equations  (I3)  emd  (l4)  will  not  be  satisfied 
Identically  by  the  vector  X  for  a  given  source  function. 

For  a  magneto- Ionic  medium  with  the  gyrotroplc  axis 
parallel  to  the  x^-axls,  the  Inverse  permittivity  matrix  will 
have  the  canonical  form. 
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^11  ^12  ° 
■^12  ^11  ^ 


(15) 


For  notation  purposes  we  define  f  =  7-^  and  use  the  subscript 
notation  to  denote  partial  derivatives  of  f. 

Combining  the  above  results  we  find  that  the  gauge  condi¬ 
tion  requires  the  following  set  of  homogeneous  equations  to 
be  satisfied  for  a  given  set  of  parameters  describing  the 
magneto-ionic  medium, 


^23  "^13  '^23 

^11 

"  o" 

^13  ^23  "^15 

^12 

= 

0 

1 

o 

CM 

CM 

r-| 

O 

_ 1 

0 

(16) 


If  non-trivial  solutions  are  to  exist  for  the  parameters 
>22*  It  Is  necessary  that  the  determinant  of  the  coef¬ 

ficient  matrix  vanish.  Therefore  we  have  the  characteristic 
equation. 


'23 

-'ij 

CM 

1 

15 

‘25 

-^5 

0 

fll+f22 

0 

Indeed  the  determinant  does  vanish  and  there  are  non-trlvlal 
solutions.  However,  it  is  also  apparent  that  it  is  necessary 
that 

(fll+f22)  ^<12 


=  0  . 


(18) 
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In  terms  of  the  vector  potential  function  this  requires 

=0  (19) 

where 


Therefore,  we  conclude  that  for  ^2.2  ^  ^  divergence  of  ^ 
must  satisfy  Laplace's  equation  in  the  plane  transverse  to 
the  gyrotropic  axis . 

However,  we  also  realize  that  if  then  the  charac¬ 

teristic  equation  is  satisfied  identically  and  the  solution  of 
(14)  represents  exact  solutions  of  Maxwell's  equations.  This 
condition,  of  course,  reduces  the  medium  to  an  isotropic  medium. 

From  the  structure  of  the  mathematical  system  and  some 
a  priori  knowledge  of  the  ionosphere  we  are  motivated  to  con¬ 
sider  the  possibility  of  obtaining  solutions  of  (14)  for  values 
of  the  plasma  parameters  which  will  give  the  conditions 


^11  ^35 


).l2-0 


(20) 


These  restrictions  are  necessary  so  that  the  imposed  gauge 
condition  yields  closed  form  solutions  representing  approxi¬ 
mate  solutions  to  Maocwell’s  equations. 


Useful  Approximations  for  the  Permittivity  Matrix  Pertaining 
to  a  Honi^eneous  Ionosphere 

Restricting  our  attention  to  the  regions  of  the  ionosphere 
and  frequency  spectrum  for  which. 
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10^  _<  f  <  10^ 

h  _>  160  km  (21) 

N  10^  electrons/c  .c . 

we  find  the  following  approximations  to  be  applicable  to  the 
elements  of  the  inverse  permittivity  matrix; 

=  -  1/X 

>ll2  =  -  i  Y/X  X  »  >>  1  >>  Z  (22) 

=  -  1/X 

X:  Normalized  Plasma  Frequency  Squared 
Y;  Normalized  Gyro  Frequency 
Zt  Normalized  Collision  Frequency 
Therefore,  we  find  that  the  conditions  stipulated  in  (20)  are 
approximately  satisfied.  This  provides  sufficient  validity 
for  the  use  of  (14)  as  determining  the  potential  function  from 
which  the  electric  and  magnetic  fields  may  be  determined  from 
(7)  amd  (8) .  It  is  noted  that  the  frequency  range  specified 
in  ( 2l)  contains  the  VLF  spectrum  so  the  approximations  in 
(21)  are  hereafter  designated  as  the  "VLF  approximations." 

The  theory  is  also  applicable  to  a  frequency  spectrum 
for  which 

f  >  10^  cps .  ( 23) 

Combining  this  restriction  with  the  altitude  range  given  in 
( 21)  we  find  that  the  elements  of  the  inverse  permittivity 
matrix  may  be  approximated  as 
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11  ~ 

^-iXY  X  «  1;  «  1;  Z  <<  1  .  (24) 

The  approximations  stated  in  (24)  shall  hereafter  be  desig¬ 
nated  as  the  "VHP  approximations." 

Solution  for  the  Vector  Potential  for  the  VLF  and  VHP  Spectrum 
The  abov$  discussion  has  provided  sufficient  validity  for 
the  use  of  the  following  set  of  equations  for  describing  the 
radiation  of  current  sources  in  the  ionosphere: 


3 

o 

(25) 

3  =  7  X  X 

(26) 

^  “  ItFe”  ^^o^  )i7(7.^)  ] 

(27) 

We  now  seek  a  solution  for  the  vector  X  in  terms  of  a 
dyadic  Green's  function.  Due  to  the  linearity  of  the  above 
equations,  we  assume  the  vector  potential  may  be  expressed 

^(x)  =J  d(x,x')?(x')dx'  (28) 

V' 

It  is  assumed  that  the  source  function  ?(x')  is  well  defined 
over  the  source  coordinates .  Combining  ( 25)  eind  ( 28)  it  fol¬ 
lows  that  the  Green's  function  must  satisfy 

(7^  +  kQK)d(x,x')  =  -  «(x-x')l  (29) 

Employing  Pourier  transform  theory  we  find  the  elements 
of  the  dyadic  Green's  function  to  be 
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-ik,R 


®11  =  87  L— R -  +  — R - J  > 


®22  "  ®11  " 


,  ^  -i  re“-^'^l^ 

’12  =  87  L-T? - 

’21  "  "  ®12  ’ 


ik,R  g-ik2R-i 

- R— J’ 


1  /e 


ik-,R 


’53  "  47' 


—  -fcrw— *% 

(^) , 


°13  "  °31  "  ®23  =  ®32  =  ^  ’ 

1  1 

ki  =  ^o^'^ll‘^^'^l2^^'  ^2  ^  ^o^'^ll*^'^l2^^'  ^^5  ~ 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


(36) 


R  =  Ix-x' 


Having  determined  the  Green's  function  for  the  vector 
potential  it  is  a  simple  matter  to  determine  the  electric  and 
magnetic  fields  from  ( 26) ,  ( 2j) ,  and  ( 28) . 


RADIATION  FROM  AN  ELECTRIC  DIPOLE  IN 
A  MAGNETO- IONIC  MEDIUM 


Radiation  Fields  of  an  Electric  Dipole  Opearting  in  the  VLF 
Spectrum 

Applying  the  VLF  approximations  to  the  above  results  with 
the  additional  constraint  that  >>  1,  the  electric  and 

magnetic  fields  are  given  in  spherical  components  as 


k^sin  O  -ik  R 

Ep  =  (-Srs - )e"'^  -  (P,  -  IP,) 


=  ( 


k^cos  ©  g-ik^^x/-Y  R 

"57“:”^®  - R - ^*^1  ■  ^*^2^ 


o 

ikjcos^©  R 

"cp  =  (“TO— -  ^"2) 


(37) 

(38) 

(39) 
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«R 


=  0 


H, 


o  /X  <.l/^^icp  e  o'  /„  ^ 

■G  ~  ‘■57“  ®  - 5 -  2^ 


cok. 


.°°=  ®  ,X  ^V2.lcp  * 

-  ^  j  e  - - - 


(40) 

i  -  ^"2) 

(41) 

.  ( Pj  -  ipj) 

(42) 

“cp  - - 5i - 

P  =  u 
n  icu  n 


u^  are  directional  cosines  specifying  the  direction  of 
the  dipole  with  the  coordinate  axes . 


VLF  Radiation  Patterns  for  an  Electric  Dipole 
Defining  the  Poynting  vector  as 

?  =  1/2  Re  X  i!*}  (45) 

the  following  components  are  obtained  for  an  electric  dipole 
operating  in  the  VLF  spectrum. 

(ldl)\^cos^e  ,Y  .1/?  .  2  2s 


®R  “ 


Sg  = 


-——7— *—7  ^  ' 

64Tr  05®  qR 


2^2 
u,  +  u- 


2  ^ 

(idl)  k-'^sin  G  cos  G  ^  t  /o  “n  “9 
'  '  o  /As  1/^  f  1 

^ ( - 2 - > 


64Tr  ooc^R 


Sep  =  0  . 


(44) 

(45) 

(46) 


Normalized  radiation  patterns  for  the  above  results  are  plotted 
in  Figure  1.  The  normalizing  factor  being 
(Idl)2(u2  +  u2)k^ 


So  = 


,9  ^ 

52ir  05® qR 


(47) 


which  is  the  maximum  value  of  the  Poynting  vector  for  the 
dipole  in  free  space. 


47 


Figure  1 — VLF  radiation  patterns  for  an  electric  dipole  trans¬ 
verse  to  the  gyrotropic  eucis. 
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Radiated  Power  for  an  Electric  Dipole  Operating  In  the  VLF 
Spectrum 

The  net  power  radiated  by  the  source  dipole  Is  defined 
by  the  Integral, 


P 


^^•n  dA 


A 


(48) 


where  A  represents  any  well  defined  closed  surface  enclosing 
the  dipole  source. 

Putting  the  results  of  (44)  through  (46)  Into  (48),  the 
net  radiated  power  for  an  electric  dipole  operating  In  the  VLF 


spectrum  Is  obtained. 

(Idl)^(ui  +  ,X  ^l/2 

^  - - TOST - 


-Y^ 


(49) 


If  we  define  the  power  radiated  by  the  same  dipole  In  free 
space  as 


(idi)  V 


**o  TTmue" 


(50) 


the  result  In  (49)  becomes. 


u,  +  u- 


(51) 


Since  the  coefficient  In  (5I)  Is  large  we  realize  that  the 
radiation  of  the  dipole  Is  enhanced  by  the  presence  of  the 

« 

Ionosphere . 


Radiation  Fields  for  an  Electric  Dipole  Operating  In  the 
f^pectrum  ^  ^ 

Applying  the  VHP  approximations  to  the  results  obtained 

In  (25)  through  (36),  Including  the  radiation  zone  restrict lor s , 

we  obtain  the  electric  emd  magnetic  fields  of  an  electric  dipole 

operating  In  the  VHP  spectrum. 
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k  sin  e  r  -i(k  R-cp)  i(]c.R+q3). 

=  -5iFc7-  (ar)L(V^P2)  IT - (V«2)  IT - ]- 

(52) 

] 


k^cos  ep  -l(k  R-cp)  g-i(k2R4CP) 


+  (P^+lPj)^ - ^ 


k^sin  e  r,«-ik,R 
47r€_ 


,  3-  p  1 
L  R  > 


-  XY(l-X+XY)sin^e]g  (Pl-iP2) 

^  [1  +  xy(  i-x-xy)  (Pi+iP2) 


2  xy  sin  e  cos  6  ^  ^  '3 

>  J  ' 


.-ik_R 


Hr  =  0  , 

-ik  05 

“©  =  "8F— 


[(1  -  X  +  xy)^/2  e 
-  (1  -  X  -  xy)^/2  e . (Pl+iP2)]  > 

k  05  cos  0  p  ^-i{k,R+cp) 

H,  •  -2-8i - [d  -  X  .  xy)V2 

.  (1  -  X  -  x^)V2  (-l-^‘’2)] 


k  05  sin  0  ^-ik,R 
o  e  3  _ 

^  R"^  ^3  • 


(53) 


(54) 

(55) 

(56) 


(57) 


VHF  Radiation  Patterns  for  an  Electric  Dipole 

The  ccmponents  of  the  Poynting  vector  are  readily  cbtair.a 
from  the  above  field  expressions. 

(Idl)  V 


S 


-  ( 


r  2  2  2  2  2  1 

o  -  V  '}  (u. +U9)cOS  0  +  u,sin^0 

L  1  2^  5  J 

( Idl) ^k^sin0cos0u 


2"-^-^:^ - cos^j^rLuj^cos(  /32R-IP)  -U2Sin(  J  , 

( 58) 


iSir  cjoe^R 
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=  - 


( Idl)  ^k^sinGu,  n  -i 

-2-^  (xy)  sin)3j^R|^Uj^sin(  i32R-«P) +U2COs02R-^)  J  > 

(59) 


IStt  odc^R 


Sjp  =  0 


(60) 


Normalized  radiation  patterns  for  the  above  results  are 
plotted  In  Figure  2.  The  noarmalizing  factor  was  given  in 
(^7)  . 


Radiated  Power  for  an  Electric  Dipole  Operating  in  the  VHF 
i^pectrum 

Evaluating  the  integral  in  (48)  for  the  fields  pertaining 


to  the  VHF  spectrum  gives  the  net  power  radiated  by  the  dipole, 
2  2 


(61) 


We  note  that  the  radiated  power  of  a  dipole  oriented  parallel 
to  the  gyrotropic  axis  is  equal  to  the  power  radiated  by  the 
same  dipole  in  free  space.  However,  if  the  dipole  is  oriented 
transverse  to  the  gyrotropic  axis,  the  power  radiated  is  one 
half  that  of  free  space. 


RADIATION  FROM  A  MAGNETIC  DIPOLE  IN 
A  MAGNETO-IONIC  MEDIUM 


Radiation  Fields  of  a  Magnetic  Dipole  Parallel  to  the  Gyro¬ 
tropic  Axis — vtJ*  Spectruun 

Omitting  the  details,  the  electric  and  magnetic  field 
components  for  a  magnetic  dipole  parallel  to  the  gyrotropic 
axis  and  operating  in  the  VLF  spectrum  are: 
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Figure  2 — VHP  radiation  patterns  for  an  electric  dipole. 
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«R 

Hep 


0, 

Mk^sin©  ^-ik  7^/-Y  R 

°  (XA)  ^ - 


2 

■BF 

2 

Mk_sin©cos© 


bir 


(xA) 


,-ik^^/3^7=Y  R 


Mk^sln^e  ,  ^-ik^AZ-Y  R 

(-rasc-xv-v)^/^  - 

o 


o 

ib-roe. 


Mk^sinGcos©  ,  ^-ik  7Jc/-'Y  R 

°  -  (X/-Y)^/2  e - - 


Mk^sine  ,  ,  _-ik„^/x/-Y  R 

^|j^(X/-Y)V2u^,eo.2e]-  “ 


(62) 

(6?) 

(64) 

(65) 

(66) 

(67) 


It  Is  noted  that  the  above  field  equations  have  the  factor 
involving  X/-Y  which  is  much  larger  than  one.  It  is  also 
interesting  to  compare  the  above  field  structure  with  that 
for  an  electric  dipole  operating  in  the  VLF  spectrum  and  trans 
verse  to  the  gyrotropic  axis . 


VLF  Radiation  Patterns  for  a  Magnetic  Dipole  Parallel  to  the 
flyrotropic  Axis 

The  components  of  the  Poynting  vector  for  a  magnetic 
dipole  peurallel  to  the  gyrotropic  axis  are  readily  obtained 
from  the  above  equations . 


2  *5  2 

M^k^sln  © 

7—2 - 2  (X/-Y) 

SUw  <J0€qR 

M^k^sin^Gcos© 

o 


^/^(cos^©  + 
(X/-Y)^/^  , 


(68) 

(69) 


s 


flp 


=  0  . 


(70) 
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Radiation  patterns  corresponding  to  these  results  are  given 
in  Figure  5.  The  normalizing  factor  being 


Mk- 


S_  = 


iS-ir  cue. 


(71) 


Radiated  Power  for  a  Magnetic  Dipole  Parallel  to  the  Gyro- 
tropic  Axis — VLF  Spectrum 

Defining  as  the  power  radiated  by  a  magnetic  dipole 
in  free  space  we  readily  find  the  power  radiated  to  be 

P  =  (Pq/4)(x/-y) 5/^(1  +  3/2y)  (72) 

Due  to  the  restrictions  of  the  values  of  the  plasma  para¬ 
meters  in  the  VLF  spectrum  we  realize  that  the  radiation 
characteristics  of  the  dipole  are  somewhat  better  than  in 
free  space . 


Radiation  Fields  for  a  Magnetic  Dipole  Parallel  to  the  Gyro- 
tropic  Axis — Vtri*  Spectrum 

For  the  VHF  spectrum  we  find  the  electric  and  magnetic 
fields  for  a  dipole  parallel  to  the  gyrotropic  axis  have  the 
following  form: 


«R 

«e 


=  0, 


Mk^slne  r  ^-ik^R-i 

—  [(i-x+jor)-^  +  (i.x-XY)_g-L.]  , 


Mk-^slnecos© 
Mk^sin^G 


[(l-X+XY)-..g-^-  -  (l-X-XY)^-g-i-J  , 
-ik^R-, 

(XX)  . 


■ik-R- 


= 


o 

IbTTUje, 


Mk^sin^G 


®G  = 


= 


IBtoc, 


[(1+XY) 


^-ik.R  ^-ik~R 


]. 


o 

tlTROC. 


(73) 

(74) 

(75) 

(76) 

(77) 


|_(1+XY  cos^G)^g-=----  (1-XY  cos  G)  ft  -■  J  , 


(78) 
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Since  the  normalized  plasma  and  gyro  frequencies  are  much  less 
them  one  it  is  apparent  that  the  coefficients  in  the  above 
expressions  may  be  simplified  somewhat;  however,  the  effect  of 
the  ionosphere  is  evident  in  the  phase  factors  of  the  propa¬ 
gating  modes . 


VHF  Radletlon  Patterns  for  a  Magnetic  Dipole  Parallel  to  the 
iSyrotr6pic  Axis 

The  Poynting  vector  for  the  VHP  spectrum  has  the  following 
components: 

M^k^sin  ©  r  ,  n 

S_  =  2 . - 7  +  cos(kj^-k2) R  +  Y  sin  20(  1-cos (kj^-k2) rJ  , 

647r  obe^R 

°  (79) 

Mk  sin^Gcos©  2 

S  - - 2— - -  (xy)^(i  +  cos(k,-k2)R)  ,  (80) 

S^TT 

Mk^sin^©  2 

S  - - 2_ - -  (XY)""  sin(k,-k2)R  .  (8l) 

647r  ooe^R 

Radiation  patterns  for  the  above  results  are  given  in  Figure  4. 

The  data  for  a  magnetic  dipole  transverse  to  the  gyro- 
tropic  2Dcis  (parallel  to  the  x^^-axis)  is  given  below  in  outline 
form. 

Radiation  Fields  of  a  Magnetic  Dipole  Perpendicular  to  the 
6yrotropic  Axis — VLp  Spectrum 


»R 

=  0  , 

(82) 

2  "5 

Mk  cos-^©  _• 

-i(k,R-qp) 

«e 

=  -§i —  (w)  1 

'  J. 

- R -  ' 

(85) 

2  2 

Mk  cos  © 

-i(k,R-?P) 

=  1  (XA) 

e  X 

- R -  ' 

(84) 

57 


=  - 


Mk^sinGcos^e  ,  --i(k,R-5P) 

-  (x/.y)1/2  =_^^1 - 


2 

iMk^cos  e 


=  - 
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VW  Radiation  Patterns  for  a  Magnetic  Dipole  Perpandlcular  to 
the  Gyrotroplc  Axis 


2  5  4 

M.  lc_cos  Q  -j,  p  o  *1 

®R - ^2 - 7  (X/-Y)^/^  [(X/-Y)  +  cos^gJ 

^  I287r  oac^R^ 
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Radiation  patterns  are  given  in  Figure  5. 
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Radiated  Power  for  a  Magnetic  Dipole  Perpendicular  to  the 
flyrotropic  Axis — Vl^  spectrian 

p  =  5/40(x/-y)5/2  p 


Radiation  Fields  of  a  Magnetic  Dipole  Perpendicular  to  the 
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Radiation  patterns  are  given  in  Figure  6. 
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BACKGROUND  STUDY  FOR  THE  CALCULATION  OF 
TRE  DRIVING -POINT  IMPEDANCE  OF  A  CYLINDRICAL 
ANTENNA  IN  OHE  ANISOTROPIC  IONOSPHERE 


Uie  derivation  of  the  driving -point  impedance  of  any 
finite  antenna  in  the  ionosphere  requires  the  calculation 
of  the  electromagnetic  fields  in  the  region  directly  adja¬ 
cent  to  the  surface  of  the  antenna.  The  quasi-static  fields 
(fields  at  a  distance  from  the  antenna  much  less  than  a 
wavelength)  must  be  determined  for  the  given  ionospheric 
conditions  surrounding  the  antenna,  and  for  the  given  shape 
and  material  of  the  antenna.  For  the  first  approach,  the 
ionosphere  is  chosen  to  be  an  infinite  homogeneous,  aniso¬ 
tropic,  lossy  (magneto -ionic)  medium.  The  anisotropy  is 
caused  by  the  earth's  magnetic  field,  whose  effect  is  assumed 
uniform  in  the  region  of  the  antenna.  The  antenna  is  chosen 
to  be  a  finitely  long,  thin,  center-driven  cylinder. 

From  a  search  of  literature  certain  references  have 
been  found  indirectly  applicable  to  the  problem  of  calculat¬ 
ing  the  quasi-static  fields  at  VLF  radiating  from  a  cylindrical 
antenna  in  the  ionosphere.  Progress  in  the  mathematical  form¬ 
ulation  of  this  analysis  has  been  made  by  Bunkln  [1];  by 
Ku^l  [2],  who  wrote  a  working  paper  based  on  Bunkin' s  analysis; 
and  by  Friedman  [3],  who  considered  propagation  from  a  dipole 
in  the  ionosphere  to  the  earth,  using  a  spherical  earth  model. 
However,  these  three  papers  considered  only  dipole  radiation 


62 


63 


far  from  the  antenna;  and  although  the  formulations  are 
instructive,  the  solutions  are  not  useful  for  practical 
applications  of  calculating  impedances  of  antennas  in  the 
ionosphere.  A  recent  mathematical  formulation  of  the  prob¬ 
lem  by  Mittra  [4]  is  more  directly  applicable.  The  best 
background  reference  is  a  book  being  published  by  McGraw- 
Hill  [5]. 

Mittra  formulated  a  mathematical  technique  by  which 
one  can  derive  the  quasi -static  Green's  function  solution 
to  the  tensor  Helmholtz  equation.  Because  of  the  ionospheric 
anisotropy.  Maxwell's  equations  are  written  as  a  set  of 
vector  equations,  which  can  be  abbreviated  in  matrix  form. 

If  Cartesian  coordinates  are  used  with  the  z-axis  aligned 
with  the  earth's  magnetic  field  and  is  suppressed, 

then  Maxwell's  equations  for  an  anisotropic  ionosphere  are 
as  follows  [6]: 

t  X  E  +  jon^H  *  5  (a) 

V  X  H  -  =  J  (b) 

e  ■ 

V-H  -  0  (c) 

7*  (eS)  -  p/Cq  (d) 

Where  S,  S.  are  the  electric  and  magnetic  field  vectors,  respec¬ 
tively; 

|i  ,  e  are  the  free-space  permeability  and  permittivity, 
respectively; 

3  is  the  source  electric  current  density; 

p  is  the  volumetric  source  charge  density; 
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0)  is  the  radian  frequency; 

e  is  the  (unitary)  relative  permittivity  matrix; 


N  is  volumetric  electron  density; 

m,  q  are  electronic  mass  and  charge  (magnitude), 
respectively; 

is  earth's  magnetic  field; 

m  are  the  plasma,  cyclotron,  and  collision 

p  q  r 

radian  frequencies,  respectively. 

The  use  of  the  above  forms  of  ej^2'  presumes  that  the 

contributions  of  the  ions  can  be  neglected  in  conqparison  to 
the  effect  of  the  electrons.  This  should  be  a  satisfactory 
assunqption  for  frequencies  of  VLF  and  above  and  for  iono¬ 
spheric  heights  above  the  E-layer. 

Substituting  (lb)  into  (la)  yields  the  wave  equation 
in  matrix  form  (special  case  of  the  inhomogeneous,  tensor 
Helmholtz  equation): 

ccf  -  c  S  =  ,  (3) 

where  C  is  the  curl  in  Cartesian  coordinates,  written  in  matrix 

2  2 

k  =  (0  |i_e_  . 
o  o  o 


form;  and 


The  form  of  the  matrix  e  in  (le)  can  be  diagonal¬ 
ized  by  a  transformation  and  consequently  will  reduce  the 
matrix  form  of  (5)  into  a  considerably  simpler  form  [4]. 

This  transformation  is  the  key  to  obtaining  an  analytical 
solution  to  (3)  in  the  quasi-static  zone  for  the  VLF  spectra. 
As  a  result  the  quasi -static  Green's  function  solutions  can 
be  derived. 

The  method  by  which  we  expect  to  derive  the  antenna 
driving -point  impedance  involves  the  following  two  steps: 

(a)  The  calculation  of  the  electromagnetic  fields 
near  the  antenna,  for  an  assvimed  current  distribution, 

(b)  Computation  of  the  driving -point  impedance 
for  an  assumed  current  distribution  using  an  "induced  emf" 
technique  or  an  iterative  (Hallen)  technique. 

At  this  time,  the  "induced  emf"  technique  appears 
simpler  and  will  be  attempted  first.  From  this  outlined 
approach,  it  will  be  necessary  to  develop  basic  theorems, 
stating  the  meaning  of  equivalence,  reciprocity,  duality,  and 
the  Poynting  vector,  in  a  magneto-ionic  medium. 

Recent  results,  which  have  yet  to  be  checked  and 
verified,  show  that  the  Green's  function  solutions  for  the 
electric  fields  are  of  the  form  of  an  operator,  0(r,z)  oper¬ 
ating  on  a  function  of  g(cj^j^,  c^^,r,z)  where  r,z  are  the  radial 
and  longitudinal  coordinates  respectively  in  a  cylindrical 

coordinate  system  whose  z-axis  is  aligned  with  the  earth's 

2  2  1/2 

magnetic  field,  and  r  =  (x  +  y  )  .  For  exan^le,  with  the 

antenna  aligned  with  the  earth's  magnetic  field  (longitudinal 
orientation),  the  z -coitponent  of  the  electric  field  is  the 
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following : 


=  -j 


J, 

az 


/  2  2  v.2\-l/2 

^®11  ®11®55^  ' 


(4) 


^is  expression  reduces  to  the  isotropic  case  by  letting 
»  1.  in  which  case  (4)  simplifies  as  follows: 


E_  = 


4T(ie, 


+  r 


2  j -1/2 


(5) 


The  equation  (5)  gives  the  correct  E  field  component  for  the 
aforementioned  antenna  located  in  free  space. 

There  is  little  theoretical  work  done  describing 
propagation  of  electromagnetic  energy  through  a  magneto -ionic 
medium  in  quasi -static  zone  in  the  VLF  spectrum.  Consequently, 
a  thorough  investigation  of  fundamental  concepts  (e.g.,  the 
aforementioned  equivalence,  etc.)  must  be  con^leted  to  insure 
that  the  analysis  is  valid.  The  study  of  these  concepts  is 
deliberate  and  necessary  for  our  deriving  valid  results. 
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AIR-TO-SUBSURFACB  VLP  PROPAGATION 


Part  I  -  Air  to  Submarine  Communications 

The  purpose  of  this  paper,  presented  in  two  parts,  is  to 
report  on  the  investigation  into  the  relative  merits  of  verti¬ 
cal  versus  horizontal  airborne  VLF  antennas  for  communication 
with  submerged  submarines  at  large  ranges.  In  this  paper  some 
qualitative  aspects  of  the  problem  are  considered,  while  the 
following  paper  is  more  quantitative.  The  physical  model  used 
here  is  a  spherical  earth  enclosed  by  a  concentric  layer  of 
air  2uid  a  sharply  bounded  anisotropic  ionosphere. 

The  basic  antenna  with  which  we  are  interested  is  a  long 
wire  trailed  behind  an  aircraft  at  some  height  above  the  earth. 
In  the  frequency  range  of  10  -  50  kc  the  radiation  pattern  of 
such  eui  euitenna  is  not  well  known;  the  wire  may  be  of  the 
order  of  a  wavelength  while  the  aircraft,  against  which  the 
wire  is  energized,  will  be  a  small  fraction  of  a  wavelength 
in  dimensions.  Therefore,  the  usual  approximation  of  a  sinus¬ 
oidal  current  distribution  may  not  be  very  accurate.  In  the 
absence  of  a  precise  radiation  pattern,  we  assume  that  the 
emtenna  may  be  replaced  by  its  equivalent  electric  dipole 
moment  which  will  in  general  have  a  vertical  as  well  as  a 
horizontal  component. 

The  field  which  the  airborne  trailing  wire  antenna  will 
cause  in  sea  water  at  large  ranges  can  be  computed  from  two 


68 


69 


different  theories — mode  theory  and  ray  theory.  These  two 
theories  are  equivalent  but  each  possesses  different  practi¬ 
cal  advantages .  To  obtain  a  gross  estimate  of  what  is  hap¬ 
pening  in  the  earth-air-lonosphere  wave  guide,  mode  theory 
is  the  easiest  for  computations  if  some  approximations  can 
be  made.  If  a  more  accurate  estimate  is  desired,  then  ray 
tracing  provides  the  capability  of  variable  ionospheric  and 
earth  reflection  coefficients,  enabling  one  to  inject  into 
the  problem  the  transition  from  day  to  night,  solar  disturb¬ 
ances,  and  nuclear  weapon  perturbations.  Both  mode  theory 
and  ray  theory  are  frequency  sensitive;  mode  theory  tending 
to  be  more  tractable  at  ELP  and  ray  theory  at  VLF.  Ray 
theory  is  used  in  the  following  discussions  since  we  are 
interested  in  a  .frequency  range  from  10  to  JO 

Some  nomenclature  that  we  use  is  illustrated  in  Figure 
1.  The  pleuie  of  incidence  is  the  plane  of  the  paper  and  is 
defined  by  the  great  circle  route  between  the  airborne  trans¬ 
mitter  emd  the  submerged  receiver.  Two  representative  rays 
between  the  transmitter  and  receiver  are  shown.  Electric 
fields  normal  to  the  plane  of  incidence  (horizontal 

anisotropic 


Figure  1 . 
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polarization)  are  denoted  by  the  subscript  "Jj "  and  electric 
fields  parallel  to  the  plane  of  incidence  (vertical  polari¬ 
zation)  are  denoted  by  the  subscript  " 

The  vertical  dipole  associated  with  the  trailing  wire 
antenna  will  radiate  only  vertically  (|  |)  polarized  electric 
fields  for  any  plane  of  Incidence.  The  horizontal  dipole 
will,  however,  radiate  electric  fields  polarized  vertically, 
horizontally,  or  a  combination  of  the  two  depending  on  the 
orientation  of  the  horizontal  dipole  with  respect  to  the 
plane  of  Incidence. 

Part  II  will  show  more  precisely  the  electric  field 
components  In  sea  water  caused  by  vertical  and  horizontal 
dipoles  at  various  heights  above  the  earth.  At  this  point, 
however,  we  merely  state  some  of  those  results  and  discuss 
•  their  possible  Implications  for  trailing  wire  antennas.  Con¬ 
sider  for  a  moment  a  receiver  In  sea  water  in  a  direction 
which  Is  broadside  to  a  trailing  wire  antenna.  The  effective 
vertical  moment  of  the  antenna  will  cause  eui  electric  field 
in  sea  water  which  Is  parallel  to  the  plane  of  Incidence. 

See  Figure  2. 


airborne 
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^  trailing  wire  Ej^ 

emtenna 

receiver  In 
sea  water 

plane  of  incidence 

Figure  2.  Top  View  of  a  Propagation  Path. 
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The  effective  horizontal  moment  will  cause  an  electric  field 
In  sea  water  which  Is  normal  to  the  plane  of  Incidence .  It 
Is  not  unreasonable  that  Ej_  and  E||  could  be  of  the  same  order 
of  magnitude,  and.  In  general,  the  phase  of  Ej^  and  B||  will  be 
different.  Thus  the  trailing  wire  antenna  will  cause  fields 
In  sea  water  which  are  elllptlcally  polarized  In  a  plane 
parallel  to  the  sea  surface.  It  therefore  seems  possible  to 
make  the  submerged  receiver  more  Immune  to  noise  by  using 
polarization  diversity.  This  phenomenon  Is  not  readily  ob* 
servable  In  the  air  above  the  sea  water  since  the  reflected 
^  field  virtually  cancels  the  Incident  Ej_  field.  However, 
the  tremsmlsslon  coefficient  for  ^  Is  2n  cos  6  which  Is  the 
same  order  of  magnitude  as  the  transmission  coefficient  for 
E||  which  Is  2n;  where  n  =  ( JoieQa'^)  and  ^  Is  the  angle  of 
Incidence  measured  from  the  vertical.  This  Indicates  that 
If  Ej^  and  E||  are  of  the  same  order  of  magnitude  then  Ej_  and 
E||  received  In  sea  water  at  large  ranges  will  be  roughly  the 
ssune  order  of  magnitude.  Ej_  will  tend  to  be  less  than  E|| 
because  of  the  factor  cos  O  In  the  transmission  coefficient 
for  Ej^,  but  this  Is  partially  offset  by  the  fact  that  the 
magnitude  of  the  reflection  coefflents  for  Ej^  at  the  alr- 
lonosphere  boundary  Is  greater  them  the  corresponding  reflec- 


It  Is,  In  general,  true  that  the  vertical  component  of 
the  trailing  wire  emtenna  will  cause  a  larger  field  strength 
In  sea  water  (at  large  ranges)  than  the  horizontal  component 
If  both  components  radiate  the  same  power.  However,  this 
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phenomenon  is  partially  offset,  as  will  be  shown  in  Part  II, 
if  the  height  of  the  trailing  wire  antenna  above  the  earth 
is  increased.  The  fields  in  sea  water  caused  by  the  vertical 
component  are  relatively  Insensitive  to  antenna  height  while 
the  fields  caused  by  the  horizontal  component  experience  a 
very  appreciable  Increase  with  antenna  height . 

It  is  also  to  be  remarked  that  fields  computed  from  ray 
theory  tend  to  receive  their  greatest  contribution  from  the 
first  few,  or  the  lowest  order,  rays.  It  is  desirable, 
therefore,  that  these  lowest  order  rays  be  incident  on  the 
earth  in  a  region  where  the  reflectivity  is  the  greatest 
possible,  e.g.,  sea  water. 

Finally,  it  is  interesting  to  speculate  on  the  possi¬ 
bility  of  increasing  the  directivity  of  an  airborne  trailing 
wire  antenna.  One  possible  method  for  doing  this  is  sug¬ 
gested  by  the  Yagi  array.  Thus,  if  a  second  aircraft 
carried  a  passive  trailing  wire  antenna,  as  shewn  in  Figure 
5,  then  it  is  possible  that  the  directivity  in  the  direction 
of  the  receiver  could  be  increased. 


airborne  A 
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main  lobe 
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receiver  in 
sea  water 


Figure  5. 


Part  II .  Electric  Field  Components  in  the  Sea 


The  method  chosen  to  obtain  the  electric  field  components 
in  the  sea  from  airborne  vertical  and  horizontal  antennas  is 
a  geometrical-optical  technique  based  on  a  derivation  by 
Bremmer  Cl].  Bremmer  employed  a  saddle  point  evaluation  in 
the  complex  plane  of  the  original  Sommerfeld-type  integral 
equations.  This  method  has  been  modified  to  include  a  re¬ 
evaluated  convergence-divergence  coefficient  valid  in  the 
neighborhood  of  the  caustic,  a  more  accurate  estimate  of 
ionospheric  reflection  coefficients,  and  a  "cut -back"  factor 
which  accounts  for  the  effect  of  diffracted  waves  for  re¬ 
ceivers  beyond  the  caustic. 

The  composition  of  the  received  electric  field  for  air¬ 
borne  transmitting  antennas  is  considered  in  the  following 
paragraphs .  Observations  of  the  resultant  field  equations 
and  results  of  the  computer  evaluation  are  discussed. 

If  a  dipole  is  considered  to  be  at  a  height,  h,  above 
an  earth  of  effective  radius,  a,  where  the  standard  atmos¬ 
phere  has  been  assumed  (a  =  4/5  times  actual  earth  radius), 
amd  ionospheric  height,  H,  two  series  of  rays  will  reach  a 
receiver  on  the  surface  of  the  earth.  One  ray  will  be  re¬ 
flected  initially  by  the  ionosphere  (unprimed  quantities) 
euid  the  other  ray  will  be  reflected  initially  by  the  earth 
(primed  quantities)  as  shown  in  Figure  1.  The  path  length 
travelled  by  the  rays  will  be  slightly  different,  or  d'  >  d 
for  h  >  0. 
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Figure  1. 


The  electric  field  at  the  receiver  can  be  considered  to  con¬ 
sist  of  a  modified  free  space  field.  In  free  space,  the  rms 
electric  field  in  millivolts  per  meter  at  a  distance,  d,  in 
kilometers  from  the  transmitter  can  be  written  as 


-ikd 


E  =  150 


mv/mAw  of  radiated  power. 


where  the  distance,  d,  is  the  total  path  length  traversed  by 

Ztt 

the  ray  and  k  is  the  propagation  constant  in  air.  The 
coefficient  ( I50)  is  obtained  by  evaluating  the  electric  field 
at  a  distance  of  one  kilometer  from  a  short  vertical  dipole 
which  radiates  one  kilowatt  over  a  flat  infinitely  conducting 
surface . 


Since,  as  shewn  in  Figure  1,  the  path  length  of  each  ray 


is  slightly  different,  the  electric  field  for  each  ray  will 


75 


Xn  general  be  of  about  the  same  relative  magnitude  but  of 
different  phase.  Since  the  phases  will  vary  with  d  and  d', 
this  will  result  In  zones  of  strong  and  weak  fields  depending 
upon  whether  the  phases  add  or  subtract.  This  analysis  can 
be  extended  to  multihop  rays,  where  the  difference  In  magni¬ 
tude  of  the  rays  Is  predominantly  affected  by  the  number  of 
reflections  by  the  Ionosphere  and  earth  since,  at  each  reflec¬ 
tion  point,  energy  Is  lost  In  transmission  across  the  boundary. 
This  meems  that  the  more  hops  a  ray  experiences,  the  smaller 
the  amplitude  at  the  dlstemt  point  compared  to  a  ray  which 
experienced  less  hops .  Eventually,  it  can  be  seen  that  a 
limit  will  be  reached  where  rays  whose  number  of  hops  are 
greater  than  soma  number  N  will  not  contribute  significantly 
to  the  total  received  electric  field  and  can  therefore  be 
neglected . 

The  coefficient,  a,  Is  defined  as  the  ratio  of  the  square 
root  of  the  cross-sectional  area  of  a  pencil  of  rays  which 
travel  linearly  from  a  tremsmltter  to  a  receiver  to  the  actual 
cross-sectional  area  of  the  rays  arriving  at  the  receiver. 

This  coefficient  takes  Into  account  the  focusing  of  the  ray 
on  reflection  from  the  Ionosphere  and  the  spreading  caused 
on  reflection  from  the  earth.  From  the  definition.  If  a 
exceeds  unity.  It  Indicates  convergence;  If  It  Is  less  than 
unity  It  shows  divergence  of  the  ray.  If  a  spherical-earth 
concentric-ionosphere  model  Is  chosen  and  the  receiver  Is 
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located  at  the  earth's  surface,  the  value  of  a  is  always 
greater  than  unity. 

Prom  the  definition  of  a,  the  convergence  coefficient 
for  the  ray  (hj)  can  be  described  by  the  following  formula 


d.  tan  T^. 
a,  =  J.  ( - ^ 

j  a  ' 


sin  © 


1/2 


As  shown  in  Figure  2,  is  the  angle  of  incidence  at  both 
the  transmitter  and  receiver.  0  is  the  angular  distance  of 
propagation  measured  at  the  earth’s  center,  dj  is  the  total 
linear  path  length  of  the  traversed  ray,  and  "a"  is  the  effec¬ 
tive  radius  of  the  earth. 


Figure  2. 


From  the 
tan 


geometry  of  Figure  2, 

(a  +  H)  sin  y r 

- 3* - 

J  (  a  +  H)  cos  ^  ® 
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Using  this  relationship,  the  convergence  coefficient  becomes 
on  evaluation 

+  H 
a  +  H 

The  convergence  coefficient  a^,  as  defined  above,  is  valid 
for  all  regions  except  in  the  vicinity  of  the  caustic,  the 
region  where  the  sky-wave  leaves  and  returns  to  the  earth  tan- 
gentially.  At  this  point,  cos  ^  ^  and  ttj  becomes  infi¬ 

nite.  This  problem  has  been  considered  by  Wait  [2]  who 
modified  the  convergence  coefficient,  valid  in  the  vicinity 
of  the  caustic  euid  beyond,  by  considering  a  third  order  ap¬ 
proximation  (Hankel  approximation)  of  the  original  Sommerf eld- 
type  integrals  rather  than  the  usual  second  order  approxi¬ 
mation  ( Debye  expansion) .  For  a  rigorous  development  the 
reader  is  referred  to  Bremmer  [l,  page  l82]  and  to  wait's 
original  paper  [2].  Application  of  this  modification  to  the 
vicinity  of  the  caustic  allows  the  geometrical-optical  cal¬ 
culation  of  the  electric  field  to  be  continuous  through  the 
caustic  zone. 

At  each  reflection  point,  on^-the  earth  or  at  the  iono¬ 
sphere,  the  ray  is  modified  in  amplitude  and  phase  dependent 
upon  cmgle  of  incidence,  the  polarization  of  the  incident 
wave  and  the  index  of  refraction.  The  ratio  of  the  electric 
field  in  the  reflected  wave  to  the  electric  field  in  the 
Incident  wave  will  be  represented  as  four  reflection  coef¬ 
ficients: 


ttj  =  ( 


a  +  H 


e 


( 


-  a  cos 


)  cos 


7T 
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The  first  subscript  denotes  whether  the  incident  electric 
field  is  parallel  (||)  or  perpendicular  (l)  to  the  plane  of 
incidence  and  the  second  subscript  denotes  whether  the  re¬ 
flected  electric  field  is  parallel  or  perpendicular  to  the 
plane  of  Incidence.  The  reflection  coefficients  of  the 
ionosphere  can  then  be  represented  as  the  elements  of  a 
2x2  matrix,  R.  The  reflection  coefficients  for  an  iso¬ 
tropic  earth,  designated  as  0(1^),  on  the  other  hand  consist 
of  only  the  terms  ||R||  and  since  |I^  and  j_R|j  are  zero  for 
an  isotropic  medium. 

Consider  the  j  ray  which  is  reflected  initially  from 
the  ionosphere  at  an  angle  equal  to  its  angle  of  incidence 
T^j.  The  j  ray  will  be  reflected  j  times  by  the  ionosphere 
and  j-1  times  by  the  earth  when  considered  at  the  earth's 
surface.  The  total  reflection  coefficient  can  be  represented 
as  a  matrix  product  of  ionospheric  reflection  coefficient 
matrices  alternating  with  earth  reflection  coefficient  matrices, 
as  given  by 


(2) 


The  total  reflection  coefficient  matrix  for  those  rays  which 
are  reflected  initially  by  the  earth  (Q_)  are  similar  in  form 
to  with  the  addition  of  a  term  d(t^j)  to  the  right  side  of 
the  product. 


(3) 
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Note  that  is  the  seune  as  except  for  slightly  different 
angles  of  incidence  and  the  term  , 

The  electric  field  that  is  of  interest  in  submarine  com¬ 
munication  is  the  horizontal  electric  field  at  the  boundary 
in  the  sea;  therefore,  if  we  pre-multiply  the  total  reflec¬ 
tion  matrix  by  T(t^),  the  transmission  coefficient  matrix, 
the  electric  field  parallel  and  perpendicular  to  the  plane 
of  incidence  will  be  determined.  The  transmission  coefficients 
as  defined  here  are  related  to  the  reflection  coefficients  by 
the  following  relationships 

T||  =  n(l  +  R||)  (4) 

and 

Tj.  =  1  +  .  (5) 

If  the  distance  between  transmitter  and  receiver  is  such 
that  there  are  rays  which  do  not  geometrically  reach  the 
receiver,  (the  disteuice  is  greater  than  the  dlsteuice  to  the 
caustic) ,  a  surface  wave  is  created  which  diffracts  along  the 
surface  of  the  earth  to  the  receiver.  The  effect  of  such  a 
diffraction  cem  be  accounted  for  by  multiplying  the  electric 
field  at  the  geometric  limit  by  a  "cut-back  factor,"  F.  This 
factor  was  developed  by  Wait  and  conda  [5]  to  account  for  the 
Influence  of  finite  conductivity  emd  earth  curvature  on  the 
field  pattern  of  receiving  antennas .  For  waves  approximately 
tangential  to  the  surface  of  the  earth,  F  degenerates  to  the 
factor  (1  +  Rg)  where  R^  can  be  Rj_  or  R||  depending  upon  the 
incident  wave  polarization. 
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The  field  patterns  of  the  transmitting  antenna,  G.(t  ,Cp), 
and  receiving  antenna,  G  (t  ,ip),  determine  the  magnitude  of 

Jm  O 

the  radiated  or  received  ray  as  a  function  of  the  incident 
and  azimuth  angles.  The  receiving  antenna  in  the  sea  is  as¬ 
sumed  to  be  a  trailing  wire  cuitenna  whose  field  pattern  is 
only  a  function  of  the  azimuthal  angle. 

At  VLF  the  ground  wave  is  either  dominant  or  of  the 
same  magnitude  as  the  sky  wave  up  to  a  distance  of  about 
1000  kilometers.  Van  der  Pol  and  Bremmer  [l]  developed 
asymptotic  expressions  for  the  electric  field  from  a  hori¬ 
zontal  and  a  vertical  dipole.  These  expressions  will  be 
used  since  they  can  be  easily  adapted  to  computer  evaluation. 
The  ground  wave  term  will  be  symbolically  represented  as 
-icp 

E  e  o 
o 

where  and  are  the  eimplitude  and  phase  of  the  ground 
wave . 

From  the  preceding  discussion,  it  is  apparent  that  the 
total  field  at  the  receiver  can  be  represented  as  the  sum  of 
a  ground  wave  and  two  series  of  sky  waves,  one  series  cor¬ 
responding  to  those  rays  initially  reflected  from  the  iono¬ 
sphere  and  the  other  series  corresponding  to  those  rays 
initially  reflected  from  the  earth.  Thus  the  horizontal 
electric  field  in  the  sea  at  the  surface  is 
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E  =  E  e 
o 


N 

+  150  Y  G^(Cp)Fa^T^Q^G^(T^,cp)^I^.^.  J. 
J=1  ^ 

+  150  2,  G^(cp)FOLjTgQgG^(T-,cp)i^T-J 
j=l  J 


(6) 


since  the  Incident  angles  at  the  ionosphere  and  the  incident 
angles  at  the  earth  are  nearly  equal  for  the  two  series  of 
rays,  for  the  hop  ray  (t^j  fit  and  trans¬ 

mission  coefficients  and  reflection  coefficients  can  be  approxi¬ 
mated  as: 

T(t^j)  «  T(t^j)  (7) 

euid 

/"  ^11 

■’'ej  )  ***  ®  '’’ej^C  0  Rj_y  * 

The  transmitter  antenna  patterns  for  the  two  series  of  rays 
are  also  related;  for  the  horizontal  dipole 

-CJ  ?)  f9) 

cind  for  the  vertical  dipole 

G^(t^,<P)  «G^(t^,cP)  .  (10) 

By  representing  the  average  path  difference  between  the  trans¬ 
mitter  amd  receiver  for  the  ray  as  Xj  and  the  average  path 
length  as  Sj  where 


•iT' 

II 

-aj) 

(11) 

J 

+  f 

(12) 
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the  electric  fields  can  be  rewritten  for  the  horizontal 
dipole  as 

N 

E.E^e-l*o+500  f;a^(*)PajTQj  a° 

J=1  J  J 

and  for  the  vertical  dipole  as 

N 

E=E^e-"'Pot500  )^G^(cp)FajTQ^ 

j=l  J  J 

The  composition  of  the  coefficients  is 

=  i  sin  kXj  +  ^  e^^^J  +  7(1  -  R|)e"^’'*j 

a.22  =  i  sin  kXj  +  ^  e^^^j  +  7(1  +  Rj_)e"^’'*J 

=  cos  kXj  +  gi  -7(1  -  R|| ) 

^22  ~  *22  * 

For  propagation  over  sea  water,  for  some  finite  altitude 
of  the  transmitter,  the  terms  (1  -  Ry)  and  (1  +  Rj_)  are 
negligible.  When  more  than  one  ray  is  considered  and  the 
transmitting  antenna  is  at  some  finite  altitude,  the  contri¬ 
bution  of  the  second  emd  third  term  of  the  coefficients  will 
be  negligible  compared  to  the  magnitude  of  sin  kXj .  The 
received  electric  field  can  then  be  taken  as 


Following  the  same  analysis  as  given  for  the  horizontal  di¬ 
pole,  the  received  electric  field  from  the  vertical  dipole  beccanes 


^  -iks 

=  E^e“^'Po+500i  j3j.(';p)FajTQG^(Tg,cp)sin  kXj  I  '  ^  J  ( I6) 

J=1  J  J 
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N 

E  =  E^e"^'*’o  + 300  ^G^(cp)FajTQG^(T^,cp)cos  kXj 
j=l 


e 

s 


Iks 


+  X 


(17) 


Discussion 

Equation  ( I6)  and  ( I7)  Indicate  that  each  ray  Is  composed 
of  a  standing  wave  which  is  a  function  of  antenna  height  and 
to  a  lesser  degree  ionospheric  height  and  propagation  distance. 
As  evident  from  the  term  sin  kXj,  the  horizontal  antenna  has 
a  definite  antenna  gain  as  the  altitude  is  increased.  For  the 
small  argument  approximation  applied  to  equations  ( I3)  euid 
(l4),  the  vertical  amtenna  is  found  to  have  the  same  height- 
gain  function  as  given  by  Bremmer  [l],  Eckersley  and 
Millington  [4],  euid  more  recently  by  Wait  the  horizontal 

amtenna  has  the  same  height-gain  function  for  the  horizontally- 
polarized  field  component  as  given  by  Bremmer. 

Considering  again  equations  (I3)  eind  (14),  as  the  trans¬ 
mitting  antenna  approaches  the  ground,  the  far  field  space 
wave  function  becomes  identical  to  the  space  wave  given  by 
Norton  [6]  emd  more  recently  by  Biggs  [7]  for  a  horizontal 
dipole  in  the  earth  near  the  surface  (except  for  the  attenua¬ 
tion  due  to  depth  of  burial) . 

It  is  interesting  to  note  that  the  conductivity  of  the 
earth  in  the  region  surrounding  the  transmitting  amtenna 
contributes  markedly  to  the  height-gain  functions  for  both 
ground-based  and  airborne  amtennas.  This  is  because  those 
series  of  rays  from  the  lower  radiation  quadramt  are  initi¬ 
ally  reflected  by  the  earth  in  the  immediate  vicinity  of  the 
antenna.  This  region  extends  out  to  approximately  a  distance 
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of  100  miles  for  an  antenna  at  30,000  feet  and  a  propagation 

distance  of  3OOO  miles .  At  VLF  and  for  small  values  of  the 

argument  (kx.  «  kh  cos  t  <  l)  in  the  height -gain  function, 

J  ® 

the  vertically-polarized  electric  field  from  the  vertical 
antenna  will  experience  a  reduction  (maximum  of  3  db)  and 
then  an  increase  in  electric  field  strength.  The  hori¬ 
zontally-polarized  electric  field  radiated  by  the  horizontal 
dipole,  on  the  other  hand,  experiences  approximately  a  linear 
increase  for  small  arguments  of  the  height-gain  function. 

The  vertically-polarized  component  radiated  by  the  horizontal 
dipole  is  dependent  upon  the  incident  angle  at  the  earth's 
surface  and  does  not  have  as  large  a  gain  as  the  horizontally- 
polarized  component.  At  a  sufficient  altitude,  the  hori¬ 
zontally-polarized  component  radiated  by  the  horizontal 
antenna  appears  to  approach  numerically  that  radiated  by  a 
ground-based  vertical  dipole. 

The  results  of  the  modified  ray  theory  for  the  spherical 
earth- anisotropic  ionosphere  model  using  the  quasl-longl- 
tudinal  approximation  to  obtain  the  ionospheric  reflection 
coefficients  (valid  for  median  latitudes)  at  a  propagation 
frequency  of  I5  kc  compare  favorably  with  available  experi¬ 
mental  data  obtained  by  Blckel,  Heritage,  and  Welsbrod  [8] 
for  ground-based  vertical  electric  dipoles  at  I6.6  kc 
measured  by  eui  airborne  receiver  up  to  distances  of  6000  km, 
see  Figure  3.  Comparison  with  actual  airborne  transmitting 
2mtennas  has  not  been  made  but  computer  results  for  the 
vertical  electric  field  component  in  air  as  shown  in  Figure  4 
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Indicate  a  definite  Increase  of  the  electric  fields  radiated 
by  the  horizontal  antenna  of  the  order  of  20  db  per  decade 
in  height  for  altitudes  greater  than  a  skin  depth  ( as  a 
measure  of  conductivity)  in  earth  in  the  vicinity  of  the 
transmitter,  while  the  fields  radiated  by  the  vertical  an¬ 
tenna,  Figure  5>  are  essentially  unchanged  with  height  for 
the  altitudes  considered  to  date,  5  miles  at  one  mile  incre¬ 
ments.  The  fields  of  the  horizontal  dipole  are  sinusoidal 
functions  of  azimuth  as  measured  from  the  longitudinal  axis 
of  the  antenna  to  the  plane  of  propagation.  The  resultant 
graphs  of  electric  field  for  the  horizontal  dipole  are  maxi¬ 
mized  with  respect  to  azimuthal  angle,  which  in  turn  deter¬ 
mines  the  aircraft  flight  path  for  a  fixed  receiver  location. 
As  the  vertical  dipole  is  omnidirectional  in  azimuth,  no 
optimization  is  required  and  direction  of  flight  would  be 
immaterial . 

The  computer  results  for  the  domineuit  electric  field 
components  in  the  sea  for  the  horizontal  (E^p)  and  the  verti¬ 
cal  dipole  (E^)  at  a  propagation  frequency  of  I5  kc  are  shown 
in  Figures  6  and  7,  respectively.  E|p  has  a  definite  gain 
with  amtenna  altitude  of  20  db  per  decade  up  to  a  height  of 
5  miles  while  E^  is  relatively  unchanged.  At  an  altitude 
of  3  miles,  a  comparison  of  E^  emd  E^  indicates  that  they 
are  of  the  same  order  of  magnitude  but  orthogonal  in  space 
with  E^  lagging  in  cunplltude  at  large  distances. 

It  C2U1  be  concluded  that,  at  VLF  emd  up  to  altitudes  of 
5  miles,  for  the  quasi-longitudinal  approximation  of  the 


in  mill i volt s/met er/kw.  radiated  E  in  mill i volt s/meter/kw.  radiated 


AIR 


Electric  Field  in  Air  (E.) 

z 

Vertical  Dipole  at  h  =  1  and 
h  =  5  miles 

f  =15  kc,  Q-L  Approximation  to 
•  the  Ionosphere 


h  =  5  miles 


Distance  in  Miles 


h  =  1  mile 


Distance  in  Miles 
Figure  5. 


in  millivolts/meter/kw.  radiated 


89 


10' 


SEA 


10“?^ 


10 


-5 


10  ■ 


'«L 

H 


10" 


Electric  Field  in  Sea  (E^-max.) 
Horizontal  Dipole  at  h  =  1  and 


■  1  .  i .  i.  I  1  III 


J - 1 _ 1— L. 


10 


10^ 


10' 


Distance  in  Miles 


Figure  6 


91 


ionospheric  reflection  coefficients,  the  airborne  vertical 
dipole  is  the  more  efficient  radiator  to  communicate  with 
submarines  (or  surface  vessels) .  However,  the  horizontal 
dipole  experiences  a  height-gain  as  the  antenna  elevation 
is  increased  which  results  in  a  horizontally-polarized 
electric  field  component  in  the  sea  which  approaches  that 
obtained  from  the  vertical  dipole. 

The  horizontal  component  of  the  airborne  horizontal 
antenna  experiences  a  definite  antenna  gain  of  approximately 
20  db  per  decade  for  altitudes  greater  than  a  skin  depth 
in  the  earth  in  the  vicinity  of  the  antenna.  The  airborne 
vertical  dipole  on  the  other  hand  is  relatively  insensitive 
to  antenna  height  at  VLF. 

At  VLF,  the  electric  field  components  tangent  to  the 
sea  surface  (the  horizontally-polarized  electric  field  and 
the  tangential  component  of  the  vertically-polarized  elec¬ 
tric  field)  are  the  only  components  of  measurable  significance 
in  the  sea.  The  total  tangential  electric  field  measured  in 
air  which  is  equivalent  to  the  tangential  electric  field  in 
the  sea  is  the  difference  in  magnitude  of  the  Incident  and 
reflected  tangential  fields.  At  VLF,  this  value  is  of  the 
order  of  2  |n|  where  n  is  the  ratio  of  the  propagation  con- 

stemt  of  air  to  that  of  sea  water  (or  earth) .  The  magnitude 

-4 

of  n  is  approximately  10  at  VLF.  The  electric  field  in 
air  normal  to  the  sea  surface  (from  the  vertically-polarized 
electric  field)  is  approximately  twice  the  incident  field  for 
all  angles  of  incidence  less  than  the  pseudo-Brewster  angle 
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(between  89°  to  90°  ^t  VLF) .  However  in  the  sea  because  of 
the  small  value  of  the  index  of  refraction,  the  normal  elec- 
trie  field  component  is  |nj  times  the  value  in  air.  In 
summary,  the  normal  component  of  the  electric  field  is  the 
only  measurable  component  in  air  near  the  earth's  surface 
and  the  tangential  component  is  the  only  measurable  component 
in  the  sea. 

An  actual  airborne  trailing-wire  antenna  is  composed  of 
both  vertical  and  horizontal  segments.  The  antenna  can, 
therefore,  be  considered  as  a  large  effective  horizontal  di¬ 
pole  euid  a  smaller  effective  vertical  dipole.  By  appropriately 
shaping  the  trailing  wire  antenna,  a  desired  ratio  of  hori¬ 
zontal  to  vertical  can  be  achieved.  As  shewn  in  Figures  6 
and  7,  the  vertical  dipole  radiates  a  larger  component  (E^) 
than  the  horizontal  dipole  (E^)  in  the  far  field.  However, 
the  horizontal  dipole  component  (E^)  approaches  E^  with  in¬ 
creased  antenna  height  and  it  is  here  suggested  that  these 
components  can  be  made  approximately  equal  for  some  receiver 
distance  resulting  in  an  eiliptically-polarized  field.  If 
by  chemce  the  phase  difference  between  the  two  orthogonal 
field  components  happened  to  be  zero,  a  relocation  of  the 
submarine  would  introduce  a  phase  angle.  The  eiliptically- 
polarized  electric  field  would  enable  a  submarine  to  orient 
the  null  of  its  floating  wire  antenna  to  best  suit  tactical 
requirements  and  also  to  use  sophisticated  receiving  techni¬ 
ques  such  as  polarization  diversity  to  obtain  better  signal 


to  noise  ratios . 
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The  field  component  (e  )  measured  in  air  near  the  surface 

z 

of  the  earth  radiated  from  an  airborne  antenna  will  in  all 
probability  be  due  to  the  vertical  segment  of  the  trailing 
wire  antenna.  Inspection  of  Figures  4  and  5  indicates  that 
if  appropriate  effective  lengths  are  assigned  to  the  hori¬ 
zontal  and  vertical  dipole,  for  exeunple  5/6  and  1/6,  res¬ 
pectively,  and  a  consteint  current  distribution  is  assumed 
over  the  length  of  the  antenna,  the  electric  field  of  the 
vertical  dipole  would  still  be  dominant.  Doubling  the 
length  of  the  vertical  segment  will  result  in  an  increase  of 
vertical  electric  field  of  6  db.  The  actual  current  distri¬ 
bution  is  quite  complex  dependent  upon  the  physical  para¬ 
meters  of  aircraft  and  antenna. 

In  conclusion,  tests  to  determine  the  electric  fields 
in  the  sea  from  a  trailing  wire  antenna  should  be  made  with 
submerged  antennas  to  determine  the  relative  magnitude  of 
Eq,  from  a  horizontal  dipole  and  E^  from  a  vertical  dipole. 

The  elllptlcity  of  the  resultant  field  should  be  verified 
and  if  it  exists  as  suggested,  should  be  further  investigated 
to  obtain  better  signal  to  noise  ratios  and  possible  inter¬ 
ference  protection. 
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INVESTIGATION  OF  PARALLEL -WIRE  AND  COAXIAL 
SIDE -LOADED  TRANSMISSION  LINES  AS  A  TRAILING  WIRE  ANTENNA 


The  purpose  of  this  study,  started  in  July,  I96I,  was 
to  determine  the  advantages,  if  any,  that  a  parallel-wire 
side-loaded  transmission  line  might  have  when  used  as  a 
trailing  wire  antenna,  and  to  suggest  a  design  for  a  practi¬ 
cal  form  of  this  type  of  antenna  for  use  in  the  VLF  spectrum. 
This  paper  summarizes  the  results  of  the  investigation. 

With  physically  long  airborne  trailing  wire  antennas 
(lengths  in  the  range  of  thousands  of  feet),  there  is  an 
obvious  mechanical  advantage  in  end-feeding  the  antenna. 

Such  an  antenna  operating  in  the  VLF  spectrum  would  probably 
be  electrically  short  (-^  <  !)•  This  would  result  in  a 
large  input  reactance  as  compared  to  the  input  resistance. 

At  high  power  levels  such  an  input  impedance  would  be  diffi¬ 
cult  to  match  to  a  transmitter.  It  can  be  shown  by  an 
approximate  method  [1]  that  for  electrically  short,  thin 
linear  antennas,  the  input  reactance  of  a  center-fed  antenna 
is  less  than  that  of  an  asymmetrically -fed  antenna  of  the 
same  length.  Thus,  the  task  of  matching  the  impedance  of  a 
center-fed  antenna  to  a  transmitter  should  be  easier  than  for 
an  end -fed  antenna. 

Theoretical  studies  by  Harrison  [2]  and  King  and 
Harrison  [5]  show  that  the  "radiation”  currents  I^j^  and  l^,- 
of  the  parallel-wire  side-loaded  transmission-line  antenna 
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of  Figure  1  have  the  same  distribution  as  the  current  on  the 
equivalent  center-fed  linear  antenna.  Thus,  this  parallel- 
wire  side-loaded  transmission-line  antenna  should  have  the 
same  radiation  properties  as  a  center -fed  linear  antenna. 

King  and  Harrison  [3]  point  out  that  such  a  structure  may  be 
of  little  practical  value  since  it  has  no  point  of  symmetry. 

In  many  situations  it  would  be  impossible  to  approximate  the 
antenna  input  voltage,  V^,  by  coupling  to  a  tank  coil  or  by 
connecting  a  transmission  line  across  the  input  terminals 
without  altering  the  circuit  appreciably.  In  the  VLF  trailing- 
wire  application,  however,  a  relatively  short  aircraft  would 
be  towing  an  antenna  several  thousand  feet  in  length,  and  thus 
it  becomes  reasonable  to  approximate  the  aircraft  and  the 
transmitter  by  V^,  a  point  discontinuity  in  the  scalar  poten¬ 
tial  field.  As  a  consequence  of  the  similarity  of  the 
"radiation"  current  distribution  in  these  two  antennas  the 
side-loaded  transmission-line  antenna  should  realize  the 
electrical  advantages  of  a  center-fed  antenna,  while  at  the 
same  time  retaining  the  mechanical  advantages  of  an  end -fed 
antenna . 

In  this  study,  the  results  of  King  and  Harrison  [3]  were 
extended  in  two  directions:  (l)  the  theoretical  work  was  ex¬ 
panded  to  cover  a  general  parallel-wire  side-loaded  transmission¬ 
line  antenna  with  an  arbitrary  number  of  side  loads  located  in 
an  arbitrary  manner  along  the  line,  and  (2)  numerical  calcula¬ 
tions  were  made  on  various  side -load  configurations,  in  order 
to  select  the  most  promising  types  for  the  trailing  wire  appli¬ 


cation. 
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The  study  of  the  general  parallel-wire  side-loaded 
transmission  line  revealed  that  its  radiation  properties 
can  be  described  in  terms  of  the  radiation  from  a  set  of 
superin^osed  single-fed  asymmetrical  linear  antennas  [4]. 

This  is  illustrated  in  Figure  2.  The  theory  shows  that  the 
"radiation"  currents,  and  1^2'  opposed  to  the  trans¬ 
mission  line  currents,  ^21'  ^12'  general 

side-loaded  antenna  have  the  same  distribution  as  the  current 
on  the  multiply-fed  linear  antenna.  It  can  be  shown  that  the 
current  distribution  on  a  multiply-fed  linear  antenna  is 
equivalent  to  the  superposition  of  the  current  distributions 
on  a  set  of  single -fed  asymmetrical  linear  antennas.  As  a 
result  of  the  linearity  of  the  system,  we  can  separate  the 
multiply-fed  antenna  into  an  equivalent  set  of  single -fed 
asymmetrical  antennas.  Thus,  by  using  a  general  parallel- 
wire  side-loaded  transmission -line  antenna,  we  should  gain 
the  same  advantages  that  were  described  above  for  the  single 
side-loaded  transmission-line  antenna. 

Another  result  of  this  study  is  the  determination  of  an 
equivalent  circuit  of  the  general  parallel-wire  side-loaded 
transmission -line  antenna.  It  was  proven  [4]  that,  for  irp’-"- 
inpedance  and  efficiency  calculations,  the  antenna  shown  in. 
Figure  5a  can  be  replaced  by  the  equivalent  circuit  shown 
in  Figure  5b.  The  equivalent  circuit  has  basically  the  same 
configuration  as  the  antenna.  The  main  simplification  is 
that  the  effects  of  radiation  from  the  two-wire  transmission 
line  can  be  taken  into  account  by  Ivimped  impedances  in  parallel 


Figure  2 
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(a)  General  side-loaded  transmission-line  antenna 


Lossless  Transmission  Lines 
(b)  Equivalent  circuit 


where  =  input  impedances  of  the  equivalent  multiply-fed 
linear  2Uitenna 


where  n  =  l,  2,  3  . . .  m. 


Figure  3. 
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with  the  side  loads,  and  the  transmission  line  itself  can  be 
considered  lossless.  The  impedances  lumped  in  parallel  with 
the  side  loads  turn  out  to  be  where  are  the  input 

impedances  of  the  multiply-fed  linear  antenna  of  Figure  2  [4]. 

The  next  step  in  the  study  was  to  use  the  equivalent 
circuit  to  investigate  in  an  approximate  manner  the  whole 
family  of  parallel-wire  side-loaded  transmission  lines  to  see 
which  particular  types  appear  most  promising  for  the  trailing 
wire  antenna  application.  Numerical  calculations  were  made 
for  two  configurations  of  the  side  loads.  The  first  was  the 
single  side-loaded  parallel-wire  transmission-line  antenna 
with  an  inductive  end-load  illustrated  with  its  equivalent 
circuit  in  Figure  4.  This  antenna  is  the  same  antenna  as  that 
shown  in  Figure  1,  except  that  the  side-load  does  not  neces¬ 
sarily  have  to  be  located  at  the  center  of  the  structure. 

By  adjusting  the  value  of  the  inductive  impedance,  Z2,  and 
the  location  of  the  inductive  side-load,  Zj^,  the  impedance, 

Z^,  which  is  the  reflected  impedance  at  point  A-A'  in  the 
direction  of  Z2,  can  be  adjusted  to  have  a  large  value  at  the 
frequency,  fj^,  shown  in  Figure  4c.  Assuming  that  the  antenna 
is  electrically  short  (2TrhA  <  1),  the  input  impedance 

of  the  equivalent  asymmetrically -driven  linear  antenna,  is 
capacitive.  By  varying  the  inductance  of  Zj^,  the  impedance 
Z^j^  can  be  adjusted  to  have  a  large  value  at  frequency  f2. 

Thus  Zg,  which  is  the  sum  of  Z^j^  and  Z^,  will  have  a  large 
value  at  fj^  and  at  f2  as  illustrated  in  Figure  4c.  The  effect 
of  the  two-wire  transmission  line  between  the  generator,  V^, 
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and  is  simply  to  shift  the  frequencies  where  becomes 
large  so  that  the  resulting  input  resistance  and  reactance 
curves  have  the  form  shown  in  the  graph  of  Figure  5. 

It  may  be  noted  that  in  the  region  of  f^  the  input 
reactance,  of  the  antenna  is  nearly  zero,  and  its  varia¬ 

tion  with  frequency  in  that  region  is  relatively  small. 

In  the  same  region,  R^,  the  input  resistance  of  the  antenna, 
is  larger  than  'th®  input  resistance  of  the  equivalent 

a symmetrically -driven  linear  antenna,  and  its  variation  with 
frequency  is  also  relatively  small.  Thus,  an  operating  point, 
exists  at  which  the  antenna  input  impedance  has  little  or  no 
reactance  to  be  balanced  out.  If  the  input  resistance  can  be 
made  sufficiently  large,  impedance -matching  to  a  transmitter 
could  be  further  simplified.  Figure  6  shows  curves  of  input 
resistance  and  reactance  plotted  against  frequency  for  an 
antenna  of  the  type  just  described.  The  length  of  the  antenna 
is  2000  feet,  the  radius  of  the  conductors  is  1  inch,  and  the 
conductor  spacing  is  28  inches,  llie  side  load  is  an  0.840 
millihenry  inductance  (with  Q  *  100  at  J&  kc)  located  at  the 
center  of  the  structure.  Obsejrving  the  curves  of  Figure  6, 
we  see  that  at  65  kc  the  input  resistance  varies  relatively 

i 

slowly  as  a  function  of  frequency  and  has  the  value  of  I5  ohms, 
in  conqparison  to  I5  ohms  for  R^,  the  input  resistance  of  the 
equivalent  center-fed  linear  antenna.  At  the  same  frequency 
the  input  reactance  also  varies  slowly  as  a  function  of  fre¬ 
quency  and  has  a  value  of  -3OO  ohms.  Ihus,  from  an  impedance 
matching  standpoint,  this  frequency  would  be  a  good  operating 
point. 


Input  Reactance 


LO' 


R -■  Input  resistance  of  equivalent  asymmetrically-driven 
”  antenna 

Figure  5.  General  characteristics  of  the  input  Impedance 

for  a  single  side-loaded  transmission-line  antenna 


in  1000  - 
ohms 


-lOOOJ 
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Another  factor  to  be  considered  is  the  power  lost 
in  the  antenna  due  to  the  finite  conductivity  of  the 
wires  and  the  finite  Q  of  the  side  and  end  loads. 
theoretical  expressions  for  the  antenna  input  impedance 
and  current  distribution  were  derived  under  the  assunption 
that  the  conductivity  of  the  wires  is  infinite.  If  it  is 
desired,  these  expressions  can  be  modified  in  an  approxi¬ 
mate  manner  to  include  the  finite  conductivity  of  the 
wires  [4].  The  efficiency  of  an  antenna  is  normally  defined 
to  be  the  ratio  of  the  radiated  power  to  the  total  power 
input  to  the  antenna,  in  order  to  simplify  the  efficiency 
calculations,  it  is  assumed  that  the  power  losses  in  the 
wires  can  be  neglected  in  con^arison  with  the  radiated  power 
and  the  losses  in  the  side  and  end  loads.  The  theoretical 
antenna  efficiency  shown  in  Figure  6  is  the  ratio  of  the 
radiated  power  to  the  sum  of  the  radiated  power  plus  the 
power  lost  in  the  side  and  end  loads,  which  are  assumed  to 
have  a  Q  ■  100  at  78  kc.  At  the  chosen  operating  frequency, 
65  kc,  this  efficiency  is  8^.  The  "overall"  antenna  effici¬ 
ency,  which  must  also  include  any  network  required  to  match 
the  antenna  inqpedance  to  the  transmitter  output  impedance, 
will  be  less  than  8^  due  to  losses  in  the  matching  network. 

nie  antenna  in  this  exanple  could  also  be  operated  at 
70.2  kc,  \diere  the  input  reactance  is  zero  and  the  efficiency 
is  about  20^.  For  this  case,  the  input  resistance  is  larger 
(although  not  as  stable  with  respect  to  frequency),  and  any 
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required  matching  network  should  be  simpler  and  more 
efficient  than  the  network  required  to  match  the  same 
transmitter  to  other  VLF  antennas,  which  may  have  a  large 
reactive  part  in  the  input  impedance. 

From  the  above  modified  definition  of  antenna  efficiency, 
it  might  appear  that  a  higher  antenna  efficiency  could  be 
obtained  with  no  side  loads.  This  is  not  the  case,  since 
it  can  be  shown  that  with  no  side  loads,  a  parallel-wire 
transmission  line  with  spacing  small  compared  to  a  wave¬ 
length  does  not  radiate  appreciably  because  of  the 
cancellation  of  equal  and  opposite  currents  in  the  wires, 
nie  side  loads  create  an  unbalance  in  the  currents,  which 
results  in  radiation  from  the  structure.  Ohe  losses  in 
the  wires  are  still  present  even  if  there  are  no  side  loads. 

The  second  configuration  chosen  for  numerical  investiga¬ 
tion  is  the  double  side-loaded  parallel-wire  transmission 
line  shown  in  Figure  7.  In  this  case  the  far  end  is  short 
circuited.  A  qualitative  description  of  the  variation  of 
the  input  inpedance  with  respect  to  frequency  is  not  as 
straightforward  as  for  the  first  case  considered.  This  is 
basically  because  the  two  side  loads  operate  as  a  coupled 
system.  Hie  theoretical  expressions,  however,  predict  that 
the  input  inqpedance  becomes  large  at  two  frequencies  which 
can  be  adjusted  by  adjusting  the  size  and  location  of  the 
side  loads. 

Figure  8  shows  the  curves  of  input  resistance,  input 
reactance,  and  antenna  efficiency  calculated  for  an  antenna 


Double  Side -Loaded  Transmission-Line  Antenna 


Equivalent  Circuit 


Figure  7, 


Double  side-loaded  transmission-line  antenna 
and  the  equivalent  circuit 


Figure  8.  Theoretical  input  impedance  and  efficiency  for 
a  double  side-loaded  transmission- line  antenna 
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of  this  type.  The  length  of  the  antenna  is  8000  feet,  the 
conductors  have  a  radius  of  l/Q  inch,  the  conductor  separa¬ 
tion  is  3* 5  inches,  and  the  side  loads  are  inductances  of 
6.1  millihenrys  (with  Q  =»  100  at  18  kc).  The  location  of 
these  side  loads  is  illustrated  in  the  figure.  Again  we 
find  a  desirable  operating  point  where  the  input  resistance 
is  100  ohms,  the  input  reactance  is  essentially  zero,  and 
the  antenna  efficiency,  neglecting  the  power  loss  in  the 
wires,  is  2^.  Even  with  high  quality  inductive  side  loads 
(Q  «  100  at  18  kc),  at  this  operating  frequency,  the  effici¬ 
ency  is  low  because  of  the  small  unbalance  of  the  transmission- 
line  currents.  It  should  be  noted  that,  while  large  currents 
may  be  flowing  in  the  transmission -line  antenna,  it  is  only 
the  net  imbalanced  current  which  produces  the  radiation? 
hence  we  can  have  a  power  loss  in  the  high-Q  side  loads 
and  no  appreciable  radiation  (and  hence  low  efficiency)  if 
the  current  unbalance  is  small. 

Experimental  verification  of  the  theory  for  a  parallel- 
wire  double  side-loaded  transmission-line  antenna  is  presented 
in  Figure  9«  l^e  experimental  model  is  scaled  down  from  the 
antenna  of  Figure  8.  We  note  from  Figure  9  that  we  have  a 
desirable  operating  point  in  the  region  of  124  me.  The  high 
frequency  is  the  result  of  frequency  scaling.  This  scaling 
is  done  in  inverse  proportion  to  the  dimensional  scaling  in 
order  to  retain  the  electrical  properties  of  the  actual  VLF 


antenna . 


Figure  9.  Comparison  of  theory  and  experiment  for  a  double 
side-loaded  transmission-line  antenna 
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It  should  be  mentioned  at  this  point  that  there  are 
many  factors  to  consider  in  choosing  an  optimum  operating 
point.  For  example,  in  Figure  8  the  frequency  for  which 
=  0  is  in  a  region  in  which  the  behavior  of  is  not 
the  best  from  the  stability  point  of  view.  Moreover, 
considerably  better  antenna  efficiencies  could  be  obtained 
at  higher  frequencies,  at  the  expense  of  larger  input 
reactance.  In  any  practical  design,  all  of  the  design 
parameters  cannot  be  simultaneously  optimized,  and  compro¬ 
mises  must  be  made. 

The  parallel -wire  transmission-line  antennas  for  which 
theoretical  characteristics  are  shown  in  Figures  6  and  8 
are  not  suggested  as  designs  for  a  practical  VLF  trailing- 
wire  antenna.  The  physical  size  of  an  antenna  with  2 -inch 
diameter  conductors  and  28 -inch  spacing  prohibits  its  use 
in  the  VLF  trailing -wire  application.  An  antenna  with 
1/4-inch  diameter  conductors  and  5,5 -inch  spacing  is  some¬ 
what  more  feasible;  however,  problems  of  maintaining  the 
required  spacing  of  the  conductors  and  preventing  the  twist¬ 
ing  of  this  antenna  in  flight  are  anticipated.  The  dimensions 
for  these  antennas  were  chosen  in  order  to  facilitate  the 
conputer  calculations  of  the  antenna  characteristics. 

The  theory  developed  for  this  class  of  antennas  shows 
that  the  input  impedance  is  relatively  insensitive  to  the 
(<t/a)  and  (lA))  ratios  (here,  t  ■  antenna  length,  a  ■ 
conductor  radius,  and  b  «  conductor  spacing),  provided  that 
these  ratios  are  very  large.  Another  theoretical  result  is 
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that  the  overall  impedance  level  is  directly  proportional 
to  ln(b/a)  [4].  Consequently,  the  input  impedance  character¬ 
istics  shown  in  Figures  6  and  8  should  be  realizable  in  an 
antenna  with  smaller  conductors  and  spacing,  provided  that 
the  original  (b/a)  ratio  is  maintained.  From  theoretical 
considerations  alone,  it  should  be  possible  to  design  a 
parallel-wire  transmission-line  antenna  with  input  impedance 
characteristics  similar  to  those  of  Figures  6  and  8  in  the 
form  of  a  parallel -wire  line  embedded  in  a  dielectric, 
with  the  outside  dimension  of  the  dielectric  being  of  the 
order  of  1/2  inch  or  less.  In  practice,  however,  the  con¬ 
ductor  diameter  would  necessarily  be  small  to  maintain  the 
required  (b/a)  ratio,  thereby  increasing  the  I  R  power 
losses  in  the  wires  for  a  given  antenna  current.  For  the 
case  of  driving  this  antenna  from  a  10  kw  transmitter, 
the  antenna  current  may  be  of  the  order  of  50  anperes  or 

more.  Such  a  current  flowing  in  small  conductors  would 

2 

lower  the  antenna  efficiency.  The  1  R  power  loss  problem 
is  decreased  with  lower  transmitter  power;  however,  the 
communication  range  is  also  decreased. 

As  a  possible  solution  to  the  problem  of  size  and  aero¬ 
dynamic  instability  which  may  be  associated  with  the  parallel- 
wire  side-loaded  transmission-line  antenna,  a  coaxial  side- 
loaded  transmission -line  antenna  is  being  investigated. 

As  a  means  of  determining  the  feasibility  of  such  an  antenna, 
relative  field  intensity  measurements  were  made  for  a 
coaxial  side -loaded  antenna  and  a  coaxial  sleeve  antenna 
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shown  in  Figure  10.  The  relative  field  intensity  with  equal 
and  constant  power  input  to  both  antennas  is  shown  in 
Figure  11.  The  field  radiated  from  a  length  of  shorted 
coaxial  cable  is  also  shown  to  indicate  that  no  appreciable 
antenna  mode  exists  for  the  coaxial  configuration  with  no 
side  load  or  slot.  Much  of  the  erratic  nature  of  the 
relative  field  intensity  data  can  be  attributed  to  the 
fact  that  the  measurements  were  made  in  a  room  with  the 
shortest  dimension  equal  to  only  a  few  wavelengths.  The 
high  frequencies  are  a  result  of  frequency  scaling.  From 
Figure  11  it  can  be  seen  that  over  an  appreciable  frequency 
range,  the  coaxial  side-loaded  antenna  radiates  nearly  as 
well  as  the  sleeve  antenna  of  the  same  length.  Since  a 
theory  for  the  general  case  of  the  coaxial  side-loaded 
antenna  was  not  available  at  the  time  of  the  experiment,  no 
attempt  could  be  made  to  predict  the  behavior  of  the  relative 
field  intensity  for  the  side-loaded  antenna. 

An  approximate  analysis  of  a  coaxial  single  side-loaded 

transmission-line  antenna  was  made  in  a  manner  similar  to 

that  made  for  the  parallel-wire  antenna.  With  the  assun^tion 

that  the  outer  conductor  diameter,  d^,  is  much  larger  than 

the  inner  conductor  diameter,  d^,  it  was  found  that  the  < 

expression  for  the  input  in^edance  has  the  same  form  as  for 

the  parallel-wire  case.  The  overall  input  inpedance  level 

do 

for  the  coaxial  case  is  directly  proportional  to  In  (g— )  for 

d^  »  d^.  The  realization  of  a  practical  coeucial  trailing 

wire  antenna  with  d_  »  d.  and  small  d  would  be  difficult 

o  1  o 
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Generator  Sleeve  Antenna  -  rg58u  (Not  to  Scale) 


Side-loaded  Antenna  -  RG38u  (Not  to  Scale) 


Figure  10 

Experimental  Models  of  a  Coeoclal  Single  Side 
Loaded  Antenna  and  a  Coeoclal  Sleeve  Antenna 
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Figure  11 


Relative  Field  Intensity  for  Coaxial  Single 
Side-Loaded  Antenna  and  Coaxial  Sleeve  Antenna 
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for  10  kw  input  to  the  antenna  because  of  I  R  losses  due 
to  the  large  current  flov;in9  in  the  small  inner  conductor. 

An  analysis  of  the  coaxial  side-loaded  transmission-line 
antenna  in  the  general  case  (i.e.,  d^  is  not  necessarily 
much  less  than  d^)  will  require  a  more  general  method  since 
the  method  of  Kulterman  [4]  becomes  too  approximate  in  the 
general  case. 

The  various  mechanical  restrictions  (weight,  space, 
tensile  strength,  stability  in  flight,  etc.)  associated 
with  the  construction  and  deployment  in  flight  of  either 
the  parallel-wire  or  the  coaxial  form  of  the  side-loaded 
transmission -line  antenna  will  impose  limitations  on  the 
range  of  electrical  properties  which  are  theoretically 
possible.  Admittedly,  the  mechanical  problems  associated 
with  a  side-loaded  transmission-line  trailing  wire  antenna 
will  be  greater  than  those  for  a  single  bare  trailing  wire 
antenna  of  the  same  length;  however,  desirable  electrical 
properties  are  theoretically  obtainable  with  the  side-loaded 
antenna  (at  least  for  the  parallel-wire  case)  for  relatively 
short  lengths. 
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EXPERIMENTAL  INVESTIGATIONS 


During  the  past  year  the  experimental  program  on  the 
project  has  changed.  Several  experiments  were  completed; 
others  were  started.  There  were  six  areas  in  which  experi¬ 
mental  investigations  played  an  important  role.  These  are: 

1.  Subterranean  Communication  Experiment  (completed). 

2.  ELF  Noise  Measurements  (closed). 

3.  Two-Layer,  Sea -Air  Model  (conqpleted ) . 

4.  Coeocial  and  Toroidal  Antenna  Impedance  Measure¬ 
ments  (active). 

5.  Side-Loaded  Transmission -Line  Antenna  (active). 

6.  Three -Layer  Model  (active). 

Several  of  the  new  areas  of  investigation  involve  the  effects 
of  the  ionosphere,  including  its  anisotropic  nature;  however, 
at  the  present  time,  there  has  been  no  attempt  to  include 
the  anisotropy  of  the  ionosphere  in  any  of  the  modeling  exper¬ 
iments.  The  theoretical  analyses  of  the  ionospheric  reflec¬ 
tions  indicate  that  it  may  be  possible  to  model  the  effects 
of  the  ionosphere  for  such  radiators  as  the  trailing -wire 
antenna  by  proper  phasing  of  the  signals  on  a  two -section 
antenna,  lliere  are  problems  confronting  the  construction  of 
the  basic  three-layer  model  which  must  be  solved  before 
detezrmining  if  this  anisotropic  modeling  is  feasible. 

Some  of  the  experiments  have  met  with  great  success, 
while  others  have  presented  formidable  problems  which  are 
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only  partially  solved.  In  order  to  facilitate  a  discussion 
of  all  of  the  experiments,  each  is  individually  discussed 
in  the  following  topics. 

1.  Subterranean  Communication  Experiment.  This  task, 
now  considered  conpleted,  sought  to  determine  the  feasibility 
of  communicating  from  a  shore  station  to  a  submerged  sub¬ 
marine  utilizing  the  earth's  crust  as  the  propagating  medium 
(see  Report  EE-8J>)  •  One  full-scale  field  experiment  was 
performed  in  two  mine  shafts  located  in  the  Manzano  Mountains, 
which  provided  a  measurement  of  the  earth's  effective  con¬ 
ductivity.  Based  upon  these  measurements  and  subsequent 
calculations,  it  became  apparent  that  a  system  utilizing  the 
earth's  cz\ist  as  the  path  of  propagation  would  not  provide 
satisfactory  ranges  of  communication. 

The  theoretical  model  employed  in  the  analysis 
mentioned  above  consisted  of  a  low -conductivity  layer  situated 
between  two  other  layers  of  higher  conductivity  (Figure  1). 
Specific  dimensions  and  conductivities  were  assumed  for  the 
various  layers  for  the  purpose  of  the  analysis.  The  low- 
conductivity  layer  was  assumed  to  be  20  kilometers  thick  and 
to  possess  a  conductivity  of  4  x  10~^  mhos  per  meter.  The 

conductivity  of  the  upper  or  sedimentary  layer  was  assumed 
_2 

to  be  10  mhos  per  meter  emd  its  thickness  was  considered 

to  be  about  1  kilometer.  Ihe  lower  layer  extending  toward  the 

center  of  the  earth  was  considered  to  have  the  same  conduc- 

_2 

tivity  as  the  upper  layer;  i.e.,  10  mhos  per  meter,  and  to 


EARTH'S  CRUST 
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^  h 

0  to  1  km  Sedimentary  Layer  (4  <  a  <  10”  ) 


h  ®  1  to  30  km  Basement  Rock  Layer 

10”^  <  ^ 
e  »»  6e 


Figure  1. 

Model  used  in  determining  propagation  of  Electro 
magnetic  Waves  in  Earth's  Crust. 
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be  of  sufficient  thickness  so  that  it  could  be  considered 
infinite  in  extent.  Assuming  the  electromagnetic  wave  to 
be  propagated  in  the  low -conductivity  (center)  layer,  the 
distance  required  for  the  wave  to  be  attenuated  by  100  db 
was  calculated  for  various  frequencies  and  for  various 
modes  of  propagation  ^®ol'  ®*ol^*  results 

of  these  calculations  are  presented  in  Table  I  and  Table  II. 

The  quantities  presented  in  Table  I  consider  the  mode 
of  propagation  to  be  which  provides  the  least  attenua¬ 
tion  of  the  possible  modes,  ™ol  were 

conputed;  however,  as  evident  from  Table  II  they  experience 
greater  attenuation  than  the  Another  advantage  of 

the  mode  is  that  the  attenuation  does  not  vary  appreci¬ 

ably  with  h,  the  thickness  of  the  low -conductivity  center 
layer,  whereas  the  other  modes  are  highly  dependent  upon  h. 


Table  I 


freq. 

100 

55 

500 

50.3 

1000 

49.4 

1500 

48.3 

5000 

42.5 

TM^^  Mode 
oo 

S 

Aj^(KM)  100  db  distance  (KM) 


158 

290 

68.2 

135 

47.0 

95.4 

38.5 

80 

25.6 

60 
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Table  II  -  TE  ,  and  TM  .  Modes 
ol  ol 


freq* 

100 

860 

500 

I8l 

1000 

105 

1500 

80.5 

5000 

55 

Aj^(KM) 

s 

100  db  distance 

656 

75.6 

129 

71.0 

68.4 

66.5 

49.5 

61.5 

28.9 

52.6 

Ihese  calculations  indicate  that  useable  signals  can  be 
received  at  considerable  distances  (e.g.,  at  100  cps  a  signal 
will  travel  290  km  before  it  is  attenuated  100  db).  The  dis¬ 
tance  through  the  basement  strata  for  which  transmission  appears 
feasible,  while  appreciable,  is  an  order  of  magnitude  less  than 
desired  for  long  range  commvinication;  and  this  does  not  take 
into  account  the  additional  losses  which  would  be  encountered 
in  penetrating  the  upper  layer  to  reach  the  submarine  receiver. 
Because  of  these  disadvantages  this  method  is  considered  im¬ 
practical  for  long  range  communications. 

The  scheme  employed  to  measure  the  effective  earth  conduc¬ 
tivity  is  new  and  useful,  and  deserves  some  consideration. 

A  pulsed  transmitter  was  en^loyed  to  obtain  high  peak  dipole 
moments.  A  receiver  enqploying  a  loop  antenna,  an  oscilloscope, 
and  a  photographic  recording  camera  was  used.  The  transmsltting 
and  receiving  units  are  Illustrated  in  block-diagram  form  in 
Figure  2.  By  measuring  the  vertical  and  horizontal  conqponents 
of  the  received  signal  near  the  earth's  surface,  a  ratio  of  the 


sin(27r*2790t) 


Receiving 

Loop 

Antenna 


Figure  2.  (a)  Block  Diagram  of  ^  kc  Pulsed  Transmitter. 

(b)  Block  Diagram  of  ^  kc  Pulse  Receiver. 
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two  can  be  formed  which  allows  one  to  determine  the  effec¬ 
tive  conductivity  of  the  earth's  crust  in  the  vicinity  of 
operation.  By  using  this  ratio,  it  is  not  necessary  to 
know  the  magnitude  of  the  transmitted  signal  since  this 
tem  drops  out  of  the  calculations.  This  new  technique  may 
prove  to  be  valuable  in  measuring  the  earth's  conductivity 
for  use  in  other  problems.  A  complete  description  of  the 
system  and  the  theory  upon  which  it  is  based  is  presented 
in  the  final  report  on  the  subterranean  experiment  (EE-85) • 

2.  ELF  Noise  Measurements.  This  experimental  investi¬ 
gation,  now  inactive,  was  to  consist  of  two  phases.  The 
initial  phase  was  primarily  concerned  with  the  development 
of  instrumentation  and  data -reduction  equipment  for  use  in 
obtaining  ELF  electromagnetic  noise  measurements.  Also  in 
this  phase  some  field  measurements  were  to  be  made  and  the 
data  reduced  in  order  to  evaluate  the  overall  system.  The 
second  phase  of  the  investigation  was  to  have  concerned  the 
improvement  of  the  instrumentation  and  data -reduction  systems 
and  the  continuation  of  the  field  measurements  and  their 
evaluation  on  a  regularly  scheduled  basis.  The  first  phase 
of  the  experiment  was  completed  and  two  reports  (EE -6 9  and 
EE -70)  were  prepared  which  describe  the  results  of  this 
initial  investigation.  Continued  work  on  the  second  phase 
of  this  experiment  would  require  considerably  more  sophisti¬ 
cated  equipment,  additional  personnel,  etc.  Since  other 
agencies  which  are  better  equipped  to  pursue  this  work  are 
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concerned  with  this  problem,  the  task  has  been  set  aside  and 
it  will  not  be  resumed  unless  there  is  a  desire  for  its  con¬ 
tinuation.  A  summary  of  the  initial  phase  of  this  experiment 
is  presented  below. 

Phase  I  of  this  task  consisted  of  the  design  and  con¬ 
struction  of  electronic  instrvimentation  for  the  purpose  of 
obtaining  electric  field  strength  measurements  in  the  3  to 
30  cps  spectrum.  Ihe  completed  system  (Figure  3)  was  placed 
in  the  field  and  recordings  were  made  over  extended  periods 
during  both  the  daylight  and  night-time  hours.  Ihe  results 
of  the  experiment,  as  evaluated  with  the  data -reduction 
system  illustrated  in  Figure  4,  indicated  the  presence  of 
the  various  mode  frequencies  in  the  5-^0  cps  region  as 
predicted  by  Schumann  and  Wait^,  although  none  were  extremely 
pronounced.  A  mode  frequency  of  about  13*3  cps  was  identified; 
however,  the  other  modes  that  should  correspond  to  the  iono¬ 
spheric  situation  resulting  in  the  13*3  cps  mode  frequency 
(8.47,  22.37,  and  29*7  cps)»  were  present  but  not  as  discern¬ 
ible  as  was  the  13*3  cps  signal.  Ihe  field  intensity  at  13«3 
cps  was  less  than  at  the  other  mode  frequencies.  A  typical 
plot  of  the  data  is  shown  in  Figure  3>  A  better  data -reduction 
system  (a  larger  number  of  narrow-band  filters)  might  have  been 


Sfait,  J.  R.,  "Terrestial  Propagation  of  Very-Low-Frequency 

Radio  Waves,"  Journal  of  Research,  NBS,  Vol.  64d,  No.  2, 
March -April,  19b0,  and 

Schiunann,  W.  0.,  "(fber  die  D^pfung  der  elektromagnetischen 
Eigenschwingungen  des  Systems  Erde -Luft -Ionosphere, " 

Z.  Naturforsch  7A>  250,  (1932). 
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Vertical 


Figure  J).  Block  Diagram  of  ELF  Noise  Receiver. 


Figure  4.  Block  Diagram  of  the  ELF  Noise  Data  Reduction 
System. 


Filter  Channel 


Average  Electric  Field  Intensity,  in  millivo Its/meter. 
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Figure  5. 


Average  value  of  observed  field  strength  during 
1000-1100  (mst)  for  period  51  July  to  8  September, 
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able  to  give  better  resolution  of  the  data.  The  equipment 
used  in  performing  the  experiment  is  being  stored  in  case 
this  task  is  to  be  re -activated. 

5*  Two -Layer  Model  (Sea -Atmosphere ) .  The  two -layer 
model  described  in  Report  EE -59#  "A  Modeling  Technique 
for  Experimental  Verification  of  Dipole  Radiation  in  a 
Conducting  Half-Space,"  is  a  successful  model  for  use  in 
evaluating  the  field  magnitude  and  field  patterns  of  various 
submerged  antenna  configurations.*  It,  of  course,  provides 
a  means  for  making  in^edance  measurements  on  the  various 
submerged  antenna  configurations;  however,  a  large  salt  water 
tank  is  available  and  is  normally  used  for  making  these 
measurements.  The  model,  consisting  of  a  wading  pool  (12  ft., 
diameter  x  5  feet,  depth)  filled  with  salt  water,  was  equipped 
with  a  set  of  frames  constructed  from  dielectric  or  insulating 
material  which  served  as  supports  for  positioning  the  trans¬ 
mitting  and  receiving  antennas  at  various  locations  and  depths 
throughout  the  pool.  The  use  of  this  configuration  allowed 
frequencies  in  the  order  of  50  mc/s  to  250  mc/s  to  be  enqployed 
in  the  model  to  represent  actual  frequencies  of  5  kc/s  to 
25  kc/s  and  model  antennas  of  10  cm  to  20  cm  length  to  repre¬ 
sent  actual  antenna  lengths  of  20  meters  to  40  meters.  The 
model,  which  was  operated  atop  a  dormitory  to  eliminate  reflec¬ 
tion  problems,  was  dismantled  after  it  had  been  used  to  make 


A  paper  covering  this  investigation,  entitled  "Experimental 
Verification  of  Dipole  Radiation  in  a  Conducting  Half  Space, " 
by  W.  E.  Blair,  will  be  published  soon  in  PGAP  Transactions. 


all  the  desired  measurements.  The  model  is  readily  avail¬ 
able  should  other  problems  require  its  use.  The  emphasis 
now  has  been  shifted  to  the  problem  of  perfecting  a  model 
which  will  provide  some  properties  of  the  ionosphere  (see 
Three-Layer  Model). 

4.  Coaotial  and  Toroidal  Antenna  Impedance  Measurements. 

The  experimental  and  theoretical  work  on  the  coaxial  and 
the  toroidal  antennas  is  still  active.  After  a  summary  of 
the  work  was  prepared  covering  these  experimental  measure¬ 
ments  and  the  corresponding  theoretical  analyses,  it  was 
apparent  that  some  of  the  assunqptions  which  were  made  in  the 
analysis  of  the  toroidal  antenna  were  not  demonstrated  to  be 
valid.  Before  the  work  is  continued  these  assumptions  are 
being  checked  out  experimentally.  A  description  of  the  two 
antennas  is  presented  below. 

The  coaxial  antenna  consists  of  a  length  of  submerged 
coaxial  cable  whose  jacket,  sheath,  and  dielectric  has  been 
removed,  allowing  each  part  to  be  exposed  to  the  water  as 

illustrated  in  Figure  6.  The  center  conductor  is  allowed  to 

^w 

extend  beyond  the  end  of  the  dielectric  at  least  ,  where 

is  the  wavelength  in  sea  water.  The  original  analysis  by 

2 

Moore  used  assumptions  that  were  equivalent  to  placing  parallel 
metallic  plates  perpendicular  to  the  cable  at  both  ends  of  the 
antenna,  with  the  plates  making  contact  to  the  center  conductor 

Moore,  Richard  K.,  "The  Theory  of  Radio  Communications 
Between  Submerged  Submarines,"  Sc.  D.  Thesis  at 
Cornell  University,  1951* 
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and  the  sheath,  respectively.  In  an  experiment  which  employed 
the  plates,  the  theoretical  and  experimental  results  seemed  in 
apparent  agreement  although  the  experimental  data  was  not  con¬ 
sistent.  A  new  set  of  measurements  is  being  made  to  see  if  the 
impedance  of  the  actual  antenna  configuration  can  be  predicted 
by  using  a  reflection  coefficient  for  the  length  of  bare  center 
conductor  that  contacts  the  sea  water.  If  this  can  be  done,  it 
will  be  superior  to  the  present  approximate  solution. 

fhe  toroidal  antenna,  initially  described  by  Anderson, 
consists  of  a  multi -turn  winding  on  a  magnetic  toroid  which  girds 
the  center  of  an  elongated  conductor.  Essentially  the  antenna 
appears  as  a  submerged  transformer  with  one  single -turn  winding 
and  another  multi -turn  winding  (Figure  7)*  Anderson^  analyzed 
the  toroidal  antenna  as  a  receiving  antenna  and  demonstrated 
that  it  possessed  some  very  desirable  characteristics.  The 

analysis  as  a  receiving  antenna  is  based  upon  an  analogous  solu- 

4 

tion  of  an  electrostatic  problem  by  Stratton  .  The  results  show 
that  the  presence  of  the  elongated  conductor  placed  through  the 
center  of  the  toroid  inqproves  the  receiving  characteristics  by 
as  much  as  ^4  db,  depending  upon  the  length  of  the  elongated 
conductor,  ^ese  representative  numbers  are  based  upon  the  open 
circuit  voltage  of  the  multi -turn  winding  and  do  not  involve  the 
internal  in^dance  of  the  antenna. 

^Anderson,  Wallace  L.,  "Submerged  Antennas,"  Technical  Memoran¬ 
dum  #2,  Engineering  Experiment  Station,  University  of 
New  Mexico,  June  12,  I96I. 

4  1 

■Stratton,  Julius  Adams,  Electromagnetic  Theory,  p.  214,  McGraw- 

Hill,  New  York,  194X1 


Magnetic  toroid 


Z; 


Conducting 

Prolate 

Spheroid 


Multi-turn 
winding  on 
magnetic  toroid 


Figure  7.  Toroidal  Antenna. 
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When  the  toroidal  antenna,  threaded  by  the  elongated 
conductor,  is  evaluated  as. a  transmitting  antenna,  the 
results  of  the  solution  of  the  electrostatic  problem  by 
Stratton  are  no  longer  applicable.  The  determination  of 
the  inpedance  is,  therefore,  a  difficult  problem  which  has 
been  attacked  by  finding  the  upper  and  lower  bounds  on  the 
impedance  without  considering  the  presence  of  the  elongated 
conductor  through  the  toroidal  axis.  As  yet,  no  satisfactory 
solution  to  the  problem,  including  the  elongated  conductor, 
has  been  found.  An  experimental  measurement  did  verify  the 
gain  relationship  existing  between  the  toroidal  antenna 
operating  with  and  without  the  elongated  conductor's 
presence.  The  publication  of  the  report  on  these  antenna 
measurements  will  be  postponed  until  these  problems  are 
solved . 

5.  Side-Loaded  Transmission -Line  Antenna.  Experimental 
work  is  being  continued  on  the  side-loaded  transmission-line 
antenna  in  an  attenpt  to  find  a  geometry  which  will  lend 
itself  to  being  towed  by  an  aircraft.  The  antenna  configur¬ 
ation  which  was  first  evaluated  was  that  of  a  parallel-wire 
transmission  line  (Figure  8).  The  large  physical  size  and 
poor  aerodynamic  properties  of  such  an  antenna  make  it  vindesir- 
able  if  not  impractical  for  conventional  use.  The  configuration 
now  being  investigated  is  cocucial  in  nature.  It  is  being 
experimentally  compared  with  a  sleeve  antenna  of  identical 
size  (Figure  9)  to  see  if  it  has  sufficient  merit  to  warrant 


Transmitter 
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a  theoretical  evaluation.  If  the  results  of  the  experimental 
test  show  that  It  does  possess  characteristics  superior  to 
the  sleeve  antenna,  the  theory  will  be  developed  and  a  practi¬ 
cal  design  attend  ted. 

6.  Three-Layer  Model  (Sea^Mr-Ionosphere) .  The  design 
of  a  three-layer  model  is  the  next  logical  step  following 
the  two-layer  model  described  above.  The  results  obtained 
from  the  two-layer  model  are  valid  only  for  a  flat  sea  when 
the  ionosphere  can  be  neglected. 

The  first  approach  to  the  design  of  a  three -layer  model 
is  the  flat-earth  configuration.  Initially  this  configura¬ 
tion  was  to  consist  of  two  layers  of  salt  water  separated  by 
a  thin  layer  of  low-loss  dielectric.  With  this  arrangement 
one  layer  of  water  represents  the  sea,  the  other  the  ionosphere; 
and  the  sandwiched  layer  of  dielectric  represents  the  air. 
Because  of  the  low  dielectric  constant  of  the  styrafoam  used 
to  represent  the  air,  the  model  is  not  large  in  wavelengths 
even  though  it  was  to  be  physically  large  for  a  laboratory 
model  (16  feet  in  diameter).  Consequently,  the  reflections 
from  the  edges  of  the  model  were  too  large  to  be  neglected. 
After  numerous  atten^ts  to  alleviate  the  reflection  problems 
on  smaller  models  were  unsuccessful,  the  configuration  was 
discarded  in  favor  of  a  different  scheme. 

If  the  dielectric  to  represent  the  air  possessed  a  large 
relative  dielectric  constant,  the  model  could  be  many  wave¬ 
lengths  in  extent  for  the  same  physical  size  mentioned  above. 
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With  this  idea  in  mind,  distilled  water  was  selected  as  the 
dielectric  to  represent  the  air.  Assuming  the  relative 
dielectric  constant  of  distilled  water  to  be  80  and  other 
parameters  the  same  of  those  of  the  air  (i.e.,  the  losses 
in  distilled  water  can  be  considered  negligible),  the  elec¬ 
trical  size  of  the  model  will  be  approximately  nine  times 
as  large  as  the  model  which  enploys  styrafoam. 

Preliminary  tests  were  performed  on  a  small  model 
consisting  of  parallel  aluminum  plates  separated  by 
distilled  water  contained  in  a  plastic  liner  (Figure  10). 

The  results  indicate  that  edge  reflections  are  not  a 
serious  problem;  however,  the  evaluation  of  these  initial 
tests  is  not  conplete.  If  the  edge  reflections  are 
negligible,  the  next  step  will  be  to  replace  the  upper 
aluminum  plate  with  a  salt  water  layer  contained  in  another 
plastic  liner.  This  configuration  may  allow  the  trailing- 
wire  antenna  to  be  modeled  and  to  take  into  consideration 
the  anisotropy  of  the  ionosphere  by  the  proper  phasing  of 
two  perpendicular  antennas. 

A  large  model  based  upon  the  use  of  water  in  each  of 
the  three  layers  can  allow  the  curvature  of  the  earth  to  be 
considered.  A  sketch  of  such  a  model  is  illustrated  in 
Figure  11.  Two  concentric  hemispherical  plastic  liners  are 
forced  to  assume  the  proper  shape  by  a  slight  hydrostatic 
pressure  produced  by  filling  the  respective  layers  to  slightly 
different  levels. 
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Figure  10.  A  cutaway  view  of  a  small  model  which  employs  distilled 
water  as  a  dielectric. 
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Figure  11.  A  possible  scheme  for  obtaining  a  spherical  configura¬ 
tion  in  the  three- layer  model. 


PLANS  FOR  1963 


During  I963  ve  propose  to  continue  in  general  the 
tasks  that  were  under  way  at  the  end  of  1962,  but  with 
particular  emphasis  on  the  ELF  portion  of  the  spectrum. 

I .  ELF/Vlf  Propagation 

In  the  area  of  ELF/VLF  propagation  we  propose  to 
continue  working  on  the  extension  of  the  theory  of  com¬ 
munications  between  stations  in  the  ionosphere,  in  the  air 
and  on  the  surface  of  the  earth  and  submerged  submarines. 

We  will  continue  to  investigate  the  effects  which  could 
enhance  ELF/VLF  propagation  and  could  lead  to  in5>rovements 
in  submarine  communication  or  to  a  better  understanding  of 
its  problems. 

Air -to -Submarine  Propagation.  Work  will  continue  on 
the  evaluation  of  electromagnetic  fields  at  VLF  caused  by 
airborne  antennas  and  received  by  submerged  dipoles.  This 
will  be  extended  to  the  ELF  portion  of  the  spectrum.  The 
effects  caused  by  an  anisotropic  ionosphere  will  be  included. 
The  main  objective  of  this  task  will  be  to  provide  a  theory 
of  computational  procedures  for  predicting  field  strengths 
in  the  sea  for  an  arbitrary  propagation  path,  arbitrary 
anteiuia  height  and  orientation,  varying  characteristics  of 
the  ionosphere  and  for  any  other  significant  parameters. 
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flie  Effect  of  the  Sea  State  on  Air -to -Sea  and  Sea- 
to-Air  Propagation.  Now  that  it  has  been  determined  that 
only  the  gross  features  of  the  sea  surface  are  required 
to  describe  such  a  surface  for  the  purpose  of  solving 
this  problem,  the  effort  will  center  on  a  search  for 
mathematical  methods  to  evaluate  the  expressions  resulting 
from  Maxwell's  equations  near  a  rough  air-sea  interface. 

At  present  the  finite  difference  method  appears  to  be  the 
best  mathematical  approach  available  for  the  solution  of 
such  expressions.  The  single  two-dimensional  case  on  which 
this  method  was  tried  checked  with  the  results  obtained  by 
others.  The  use  of  this  method  will  be  extended  to  rougher 
surfaces  and  to  three-dimensional  problems,  although  it  is 
not  yet  at  all  clear  as  to  whether  this  can  be  accomplished 
or  how  it  could  be  done.  Such  calculations  are  expected  to 
be  quite  involved,  and  extensive  use  of  computers  is  antici¬ 
pated. 

Submarine -to -Air  Propagation.  In  the  initial  study  of 
the  sea -air  ionosphere  problem  at  ELF,  ideal  conditions  were 
assvuned  (flat  boundaries,  isotropic  ionosphere,  etc.),  and 
the  field  strengths  in  the  sea,  resulting  from  submerged 
electric  vertical  and  horizontal  transmitting  dipoles,  were 
confuted.  The  results  obtained  were  restricted  to  very  large 
ranges  to  satisfy  the  constraints  in  approximating  some  of 
the  integrals.  It  appeared,  at  that  time,  that  some  kind  of 
a  numerical  process  could  be  used  to  remove  these  limitations. 
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An  attempt  was  made  to  determine  the  feasibility  of  solving 
numerically  the  integrals  associated  with  a  submerged  electric 
dipole.  The  preliminary  investigation  showed  that  the  prob¬ 
lem  was  much  more  con^licated  than  at  first  anticipated. 

It  has  been  shown  that,  although  these  integrals  can  be 
evaluated  numerically,  a  considerable  effort  would  be  re¬ 
quired.  It  appears  that  other  approaches  to  the  problem  are 
available  which  should  be  superior  and  more  direct  than  the 
evaluation  of  the  integrals.  It  is  proposed  to  investigate 
these  approaches  further. 

Satellite -to-Submarine  Propagation.  It  is  assumed  that 
a  satellite  may  be  used  to  propagate  ELF/Vlf  to  the  earth 
either  from  the  ionosphere  or  from  beyond  the  F-layer  of  the 
ionosphere.  To  investigate  the  feasibility  of  this  method 
of  communication  it  is  proposed  to  continue  the  study  of 
propagation  from  or  beyond  the  ionosphere  to  points  on  the 
earth's  surface  and  into  the  sea.  To  accomplish  this,  the 
recently  conpleted  work  for  obtaining  the  far-zone  VLF  radi¬ 
ated  power  for  current  sources  in  a  homogeneous  anisotropic 
ionosphere  will  be  extended  to  the  near  field  and  to  the 
air -ionosphere  interface.  Particular  emphasis  will  be  placed 
on  the  extension  of  this  theory  to  ELF. 

II .  ELF/Vlf  Antennas 

In  the  area  of  antenna  studies  we  will  continue  looking 
for  new  approaches  to  the  development  of  antennas  for  communi¬ 
cation  with  or  between  submarines.  Particular  attention  will 
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be  devoted  to  those  antennas  that  offer  the  most  promise 
in  the  enhancement  of  submarine  communication  in  the  ELF 
portion  of  the  spectrum.  Some  attention  will  be  given 
to  the  establishment  of  controlled  directivity  of  totally 
submerged  antennas. 

Submerged  Antennas.  The  theoretical  investigation 
of  fundamental  properties  and  limitations  of  submerged 
electric  and  magnetic  sensors  will  continue  and  will  be 
extended  to  electric  and  magnetic  antennas  as  receiving 
elements.  It  is  anticipated  that  this  task  should  deter¬ 
mine  whether  directivity  can  be  obtained  with  a  submerged 
antenna  or  prove  that,  on  the  basis  of  the  current  level 
of  the  state  of  the  art,  submerged  directivity  is  not 
practical. 

Air -to -Submarine  Antennas.  The  effort  on  this  task 
will  continue  on  the  investigation  of  side -loaded  trans¬ 
mission  lines  used  as  trailing -wire  antennas.  It  is 
proposed  to  extend  the  experimental  investigation  of  the 
coeucial  form  of  the  side-loaded  transmission  line  in  order 
to  establish  a  better  correlation  between  the  radiating 
properties  and  the  input  impedance.  Previous  attenpts  to 
determine  this  correlation  have  been  hankered  by  masking 
effects  of  reflections  from  walls  and  an  undetermined  amount 
of  radiation  from  the  model  due  to  an  unloaded  side-slot. 

It  should  be  possible  to  obtain  the  desired  correlation  by 
means  of  further  input  in^edance  and  field  intensity  measure¬ 
ments  at  a  higher  power  level  in  an  outdoor  open  area,  it 
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should  he  possible  to  determine  from  these  measurements 
whether  a  more  detailed  analysis  of  the  coeucial  side- 
loaded  transmission-line  antenna  is  justified.  Since 
experimental  work  using  a  single  bare  wire  as  an  airborne 
trailing -wire  antenna  is  now  under  way,  it  is  proposed 
that  a  theoretical  study  of  the  current  distribution, 
input  impedance  and  efficiency  of  such  an  antenna  be 
undertaken  and  the  results  compared  with  those  obtained 
on  transmission -line  antennas.  An  investigation  of  the 
use  of  a  second  trailing -wire  antenna  as  a  director  should 
also  be  undertaken. 

ELF  Antenna  for  Launching  Waves  in  Air.  The  search 
for  a  suitable  antenna  for  launching  ELF  signals  in  air 
will  continue  in  order  to  find  the  best  method  for  reliably 
getting  signals  from  the  shore  to  a  submarine  at  long 
ranges.  Calculations  of  impedance,  current  distribution, 
the  effect  of  earth  conductivity  and  permittivity,  as  well 
as  field  strengths,  should  continue.  If  the  initial  study 
shows  promising  results,  limited  verification  of  the  theory 
should  follow.  It  appears  that  the  toroidal  type  of  antenna 
buried  in  a  conducting  soil  offers  promise,  and  it  is  pro¬ 
posed  to  continue  investigating  this  type  of  an  antenna. 

Antennas  in  or  beyond  the  Ionosphere.  Theoretical 
investigation  to  determine  practical  expressions  for  the 
driving-point  in^edance  of  an  electric  antenna  at  VLF  in 
the  ionosphere  will  continue.  It  is  anticipated  that  the 
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impedance  will  show  the  effects  of  the  antenna's  orienta¬ 
tion  with  respect  to  the  earth's  magnetic  field.  The 
initial  task  will  be  on  a  bare  antenna  in  a  homogeneous 
anisotropic  ionosphere.  Once  this  problem  is  solved  it 
will  be  extended  to  an  insulated  antenna  in  an  inhomogen¬ 
eous  anisotropic  ionosphere.  The  effect  that  a  sheath  of 
electrons  near  the  antenna  may  have  on  antenna  performance 
will  also  be  considered. 

III.  Experimentation,  Modeling  and  Testing 

When  all  the  variables  are  considered,  the  problem 
of  ELF/VLF  propagation  and  antenna  design  is  one  of  great 
complexity.  To  obtain  experimental  data  at  this  frequency 
is  time-consuming  and  costly.  When  conditions  such  as 
these  prevail  it  is  often  helpful  to  scale  the  experimental 
set-up  so  that  a  simulated  picture  of  the  problem  can  be 
obtained.  This  of  course  leads  to  modeling.  It  may  be 
possible  in  a  model  to  represent  to  some  extent  the  sea, 
the  air  and  the  ionosphere  and  their  effects  on  electro¬ 
magnetic  propagation  and  antenna  performance.  Even  though 
such  models  are  not  to  exact  scale  they  are  often  helpful 
to  the  understanding  of  the  problem  as  they  may  indicate  the 
order  of  magnitude  of  one  parameter  when  another  is  made  to 
change.  Often  the  direction  of  such  a  change  may  be  signifi¬ 
cant.  As  is  true  of  any  model,  the  analogy  generally  holds 
only  for  a  restricted  range  of  operating  conditions,  but 
within  such  a  range  it  may  give  useful  answers.  This  is  the 
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reason  vrhy  modeling  is  considered  so  important  in  a  study 
of  ELF/VLF  propagation  and  antennas. 

We  now  have  a  useful  technique  to  simulate  two-layer, 
sea -air  problems;  however,  extrapolating  such  experimental 
data  to  the  three -layer  problem  (sea -air -ionosphere)  has 
presented  a  number  of  difficulties.  For  this  reason  we 
propose  to  continue  working  on  the  three-layer  model, 
feeling  that  the  study  of  the  model  itself  may  help  us 
understand  some  of  the  difficulties  involved  in  ELF/Vlf 
propagation  and  antenna  investigations. 

When  the  requirement  arises  we  propose  to  model,  if 
this  is  at  all  possible,  antennas  buried  in  the  ground  or 
submerged  in  the  sea,  and  any  novel  type  of  an  antenna. 

We  feel  that  this  may  help  in  understanding  the  problems 
connected  with  submarine  communications  and  provide  indica¬ 
tions  as  to  possible  avenues  of  attacking  them. 
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